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Abstract 



^^ We study the dynamics of the five-parameter quadratic family of volume-preserving dif- 

1^^ feomorphisms of M.^ . This family is the unfolded normal form for a bifurcation of a fixed point 

\^^ with a triple-one multiplier and also is the general form of a quadratic three-dimensional map 

^^ with a quadratic inverse. Much of the nontrivial dynamics of this map occurs when its two 

r> fixed points are saddle-foci with intersecting two-dimensional stable and unstable manifolds 

f^ that bound a spherical "vortex-bubble" . We show that this occurs near a saddle-center- 

OO Neimark-Sacker (SONS) bifurcation that also creates, at least in its normal form, an elliptic 

^^ invariant circle. We develop a simple algorithm to accurately compute these elliptic invariant 

^ circles and their longitudinal and transverse rotation numbers and use it to study their bi- 

furcations, classifying them by the resonances between the rotation numbers. In particular, 
rational values of the longitudinal rotation number are shown to give rise to a string of pearls 



C^ that creates multiple copies of the original spherical structure for an iterate of the map. 

1 Introduction 

The area-preserving Henon map |Hen69] is the universal form for quadratic diffeomorphisms of 
the plane and provides perhaps the prototype for conservative systems with chaotic dynamics. As 
such it has deservedly been the subject of much study. We believe that the natural generalization 
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of this map to three-dimensions is the volume and orientation-preserving diffeomorphism 

( x + y 

y + z-e + fiy + P{x,y) 
\ z-e + ny + P{x,y) 



f{x,y,z 
where P is the quadratic form 



P{x, y) = ax + bxy + cy . (2) 

Just as the Henon map arises as the normal form for a saddle-center bifurcation, the map / 
(with P a polynomial of any degree) is the normal form for the equivalent bifurcation in M^: a 
triple-one multiplier [DM08J . The Henon map and its generalizations to polynomials of arbitrary 
degree give a basis for the group of polynomial automorphisms of the plane |FM89j . While such a 
decomposition is not known for higher dimensions, the map / is the universal form of a quadratic 
diffeomorphism with quadratic inverse [L.M98 . More generally the normal form / has an inverse 
with the same degree as the polynomial P [D MOSj . Finally, both the Henon map and / arise as 
normal forms for certain homoclinic bifurcations |GMO06] . see ^ 

In this paper we will study some of the dynamics of the map (II|. Numerical evidence 
presented in ^ will show that the map has a nonzero measure of bounded orbits primarily near 
the simultaneous saddle-center and Neimark-Sacker (SONS) bifurcation that occurs when e = 
and —4 < // < 0. We will study the normal form for this bifurcation in ^ The bounded orbits 
of / that appear in such a regime are built around the skeleton formed from its two saddle-focus 
fixed points, their stable and unstable manifolds, and the elliptic invariant circle that is created 
when this bifurcation is supercritical. 



As we recall in ^5.1, the structure of this "saddle-center-Hopf bifurcation is well-known for 
the case of a volume-preserving flow |Bro81| IHol84j . For the supercritical case, the bifurcation 
creates a vortex-bubble structure analogous to the Hill's vortex of fluid mechanics |Lam45| or 
the spheromak configuration of plasma physics | GICH80] . There are two saddle-focus equilibria 
whose two-dimensional stable and unstable manifolds coincide forming a sphere. The interior 
of the sphere is foliated by a family of two-tori enclosing an invariant circle that is normally 
elliptic. For the map, the normal form of the "saddle-center-Neimark-Sacker" bifurcation is no 
longer integrable, but it still has two saddle-foci, a Cantor family of tori, and an elliptic invariant 
circle, see §5.2[ 

The quadratic map M is approximately described by this normal form for small e and 
away from low-order resonances. We observe numerically that many of the invariant two-tori 
and the central invariant circle appear to persist for moderate values of e. In this paper we 
will concentrate on the persistence and bifurcations of the elliptic invariant circle. We develop 
an algorithm to accurately compute this circle in ^ When it is elliptic, the circle has two 
rotation numbers w = {u;l,u;t), longitudinal and transverse, respectively. We compute these 
and compare the numerical results with the normal form calculations. 

Resonances of the form m ■ lj = n lead to bifurcations that may give rise to new elliptic 
invariant circles and may result in the destruction or change of stability of the original circle. 
Interestingly, there are two types of doubling (or m-tupling) bifurcations, one in which the new 
invariant circle is a single circle that winds multiple times around the original circle (like in a 
flow), the other one in which k new invariant circles appear that are mapped to each other. The 
latter case appears when mi and m2 are not coprime. 

Perhaps the most interesting of these bifurcations we call a string of pearls. It occurs when 
iVL becomes rational and typically results in a new SCNS bifurcation for some power of /, see 



^ When this bifurcation is supercritical, the invariant circle is replaced by a set of small 
vortex-bubbles or pearls linked by nearly coincident one-dimensional invariant manifolds of a 
pair of saddle- focus periodic orbits, the string. Recently a similar bifurcation has been observed 
for dissipative 3D maps jBSV08| . 



2 Contexts 

The quadratic map (II]) arises naturally in at least three contexts. 

Suppose that / : M'^ ^ M'^ is a smooth, volume and orientation-preserving map, det Df{(^) = 
1, where S, = (x, y, z). li^* is a point on a period-n orbit, f"'{^*) = ^* , then the multipliers of this 
orbit must satisfy A1A2A3 = 1. Consequently, the characteristic polynomial of the linearization 
about the periodic orbit. 



det(A/ - DnC)) = X^-TX^ + aX-l 



(3) 



contains two parameters: the trace r and second trace a. There are eight stability regimes in 
the (r, a) plane, see Fig. m The boundaries between these regimes contain two codimension- 
two points; one corresponds to a triple multiplier {1,1,1} at (r, o") = (3,3) and the second to 
multipliers {—1,-1,1} at (t, o") = (-1,-1). These two cases form organizing centers for the 
dynamics in all eight regimes. Notice that in the interior of each of the eight regions, the fixed 
point is hyperbolic and unstable. 




l<X,=X^<o 



Figure 1 : Classification of tlie eigenvalues for a three-dimensional, volume-preserving map as a function of the 
trace r and second trace a. The eight insets are the complex planes showing the multiplier configurations relative 
to the unit circle. 

If the Jacobian, D/"(^*) has a multiplier A = 1 with algebraic multiplicity three, then A 
generically has geometric multiplicity one; consequently, the Jacobian is similar to the Jordan 
block 

/I 1 OX 

J = 1 1 . (4) 

\0 1 



We showed in |DM08j that near ^*, /" is formally conjugate to the normal form 

{x, y, zY = J{x,y,z + p{x, y, e, /i))"^ , 
p{x, y, e, fi) = -e + /iix + /X2y + a(e, fi)x'^ + b{e, fi)xy + c(e, /x)y^ + . . . , 

where the . . . indicates higher order terms in x and y only. The normal form is "formal" in 
the sense that if / is expanded in a power series to any finite degree d in the variables ^ and a 
set of sufficiently general parameters, then this degree-d map is conjugate to d5|. Generically, 
one of the two parameters (/ii,/U2) in ^ can be eliminated, see Appendix |A| leaving two 
unfolding parameters; for example if a 7^ then /ii can be set to zero, leaving the two parameters 

(e,/i = /i2). 

The normal form provides a remarkable simplification of the full map since all of the non- 
linearity is contained in the single polynomial p that depends upon only two of the variables. 
To lowest nonlinear degree, the map ^ is quadratic and if we view (a, b, c) as independent 
parameters it becomes (II). 

The quadratic map |I) also arises in the study of polynomial diffeomorphisms: it was shown 
in |LM98j , that any quadratic diffeomorphism of M^ that has a quadratic inverse and nontrivial 
dynamics is affinely conjugate to the shift-like |.BP98, map 

(x, y, z) h^ {-e + Tx + ay + z + Q{x, y),x, y) , 

where Q(x, y) is a quadratic form. This map is linearly conjugate to the quadratic case of ^ 
under the orientation reversing transformation ^ 1— > U^ where 



U 



providing that Q{x,y) = p{y,x — y), t = fi2 + 3, and o" = ^1 — /i2 — 3. Thus the map ^ is 
also a normal form in the sense of polynomial automorphisms. However, it is not known how 
to generalize this result to cubic or higher degree — unlike the planar case that was treated by 
Friedland and Milnor |FM89j . Quadratic diffeomorphisms may have inverses of degree up to 
four; these were classified for C^ in |FW98| IMaeOlj . Our map ([5| does have an inverse of the 
same degree for any polynomial P. 

The map ([5]) also appears in a third context as the normal form for certain homoclinic 
bifurcations |GMO06] . In particular consider a three-dimensional map with a saddle-focus fixed 
point such that the saddle value — the product of the multipliers of the fixed point — is 1; this, 
of course, is true for the volume-preserving case. Suppose that the one-dimensional invariant 
manifold of the saddle-focus has a quadratic tangency with its two-dimensional (spiral) manifold. 
Unfolding this singularity leads to a map of the form (M that describes the dynamics of the return 
map in a neighborhood of the homoclinic point. There are two distinct cases: when the saddle- 
focus is of type (2, 1) (or type- A |Cow73| ) . with a two-dimensional (spiral) stable manifold and 
a one- dimensional unstable manifold, then the homoclinic normal form has P{x, y) = a(x -|- y)^, 
so that the quadratic form (p| is a perfect square. When the saddle-focus is of type (1,2) (or 
type-B), then P(x,y) = ax"^. 

An additional parameter can be introduced in M to make the Jacobian deiDf = b ^ 1. 
Maps similar to this often have strange attractors and have been studied in [ACS8"3l IAA841 
l(X)STn5llDLMn6] . 

These considerations support our assertion that the normal form ^ is an appropriate gen- 
eralized Henon map for three-dimensional dynamics. 




3 Fixed Points 

We begin the study of the dynamics of ([l]) looking at its fixed points; more details are given in 
Appendix [C| 

When ea < the quadratic map (II| has no fixed points; when ea > it has two at 

e± = (x±,0,0), x± = ±,/^. (6) 

V a 

Thus a pair of fixed points is created upon crossing the saddle-center line e = 0. Stability of these 
fixed points is determined from the general stability diagram, Fig. m by computing (r, a) for the 
fixed points, see Appendix [CJ For fixed (a / 0, b, c), the stability of each fixed point depends on 
the two parameters (e, ^). When a > the stability diagram of a;+ is a diffeomorphic copy of the 
half plane t > a; the corresponding diagram for x- covers the half-plane t < a, see Fig. |2j In 
this figure, the blue curves correspond to the fixed point x+ and the red curves to x_. Note that 
over most of the parameter range, when one of the fixed points is type (2, 1) (two-dimensional 
unstable manifold), the other is type (1,2) (two-dimensional stable manifold), see Appendix [C| 
When e = and —4 < ^ < 0, the fixed points are created with a pair of multipliers 
X± = e='=-^'^*'^o on the unit circle with rotation number 



^oifJ-) = - arcsin J-- . (7) 

As we will see next, this range of ^ seems also to correspond to the range for which there are 
other nontrivial bounded orbits for the map. 

4 Bounded Orbits 

When the quadratic form P of (pi) is positive definite, it is not hard to show that all of the 
bounded orbits of (II| are contained in a cube, see Appendix IbJ To find regions of parameter 
space that have nontrivial bounded orbits, we iterated a grid of initial conditions, declaring an 
orbit to be "bounded" if it remains in a ball of radius 10 for N iterations. The resulting fraction 
of bounded orbits as a function of the parameters (e, fi) is shown in Fig. ^ For this figure, as 
for most of the numerical studies reported below we fixed the parameters of P, choosing a = 1, 
b = c=l 

Our numerical studies indicate that the domain containing bounded orbits shrinks to zero 
as ^/e. To reflect this, we considered initial conditions on a three-dimensional grid in the box 

B = [-2x+,2x+] X [-4x+,4x4.] x [-8x+,8x+] , (8) 

where x+ is the position of the fixed point ([6]). In Fig. Islwe plot the fraction of bounded orbits 
with A^ = 10 , choosing initial conditions on a cubical grid of 500^ points that covers B. It 
appears that there are bounded orbits only when e > and —4.1 ^ yu < 0, though we have not 
exhaustively searched outside this parameter domain. We also observe that, for the parameters 
of Fig. |3J there appear to be essentially no bounded orbits when e > 0.5. 

In the figure, the dark blue regions correspond to "no" bounded orbits and the darkest 
red to the maximal fraction, 4.9%. There is considerable "resonance tongue" structure in the 
figure. Some of this is similiar to Fig. [2] for the multipliers of the fixed point. In particular the 
bounded fraction appears to be nearly zero at the doubling and tripling points (e, fi) = (0, —4) 
and (0, —3), respectively, and is small near the quadrupling point (e, fj,) = (0, —2). Indeed if we 




Figure 2: Stability diagram for the map ([T| with a = 1, b = 0.5 and c arbitrary. The solid curves correspond to 
a double eigenvalue, the dashed curves to period doubling, and the dotted curves to a complex conjugate pair at 
re with lo as indicated. Bifurcation curves for the fixed point x+ are blue and those for X- are red. 



increase N to 10^ and double the grid resolution then we found no bounded orbits in B whenever 
e > 0.001 for /i = -3 and at most 0.043% when fi = -4. 

Several constant-^ slices through the full dataset are shown in Fig. |4] Note that as e ^ 0"*" 
the bounded fraction appears to converge to a smooth curve, shown as the dashed line, except 
for the resonances near ^ = 0, —2, and —3. This curve will be derived in ^5.3 



Much of the structure of Fig. [3] and Fig. [4] can be attributed to the orbits trapped by the 
two-dimensional manifolds of the fixed points, as we discuss below. 
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Figure 3: Fraction of bounded orbits in B for the map (fTl) with nonhnearity (pi and a = 1.0 and b = c = |. 
Orbits are "bounded" if they stay inside the ball of radius 10 for A^ = IC* iterates. For this calculation, the map 
was iterated a total of 1.69 x 10^^ times over 500'^ x 500^ initial conditions. The maximum number of trapped 
orbits was 6139040, corresponding to 4.91% of the initial conditions in B when (e, /i) = (0.00078, —2.73). When 
e — 0.4 there are no bounded orbits except near /i — —2.26 and —3.95 where less than 2 x 10"''' of the initial 
conditions are bounded. 



5 Saddle-Center-Neimark-Sacker Bifurcation 

In most of the regions with bounded orbits in Fig. p] the fixed points ([6]) are spiral foci with two- 
dimensional stable and unstable manifolds, respectively. These manifolds appear to intersect 
transversely for many parameter regimes and, roughly speaking, enclose a ball that appears to 
contain all of the bounded orbits |LM 98] . The structure is analogous to a vortex bubble in an 
incompressible fluid |iLam45l IMac94] . Such a flow is generic for a volume- preserving vector field 



near a saddle-center-Hopf bifurcation, as we review in ^5.1 The topology of the intersections 
changes as the parameters vary, but most often seems to include infinite spiral curves that 
connect the fixed points. This is necessarily true for a flow since every intersection point lies on 
a homoclinic orbit; however, for a map a number of different topological types of intersections 
can also occur |LMOO| . 

The most prominent bounded orbits correspond to a family of invariant two-tori enclosing 
an invariant circle that lies approximately in the plane x = 0, see Fig. [5J If the rotation number 
([7]) of the fixed points is not close to a low-order rational number, then as e tends to zero this 
structure limits on a ball that shrinks to zero as ^/e. The stable and unstable manifolds in this 




Figure 4: Percentage of the box B that is occupied by bounded orbits as a function of /i for several values of e. 
The dashed curve is the theoretical result (23 1 that is valid as e -^ 0^, away from the major resonances. 



limit appear to nearly coincide and the family of invariant tori nearly completely fills the interior 
of the ball: the structure appears to be nearly "integrable." 




Figure 5: Cutaway view of several orbits of the map (fT| for a = 1.0, 6 = c = 0.5 and (e, /j.) = (0.01, —2.4). Also 
shown are the two-dimensional stable (blue) and unstable (red) manifolds of the fixed points 

Since the bounded orbits of the map (IT]) occur primarily along the saddle-center-Hopf line, 
we discuss in §5.2| the normal form for a map near such a bifurcation. We will show that both 
supercritical and subcritical bifurcations can occur. In the supercritical case, the normal form 
creates a pair of fixed points whose two-dimensional stable and unstable manifolds intersect and 



enclose a sphere that generically contains a Cantor-family of invariant two-tori. In ^5.3 we will 
transform ([I]) into normal form, thereby obtaining the relationship between its parameters and 
those of the normal form. This will demonstrate that this bifmxation occurs on the saddle- 
center-Hopf segment of Fig. |2} 

Before proceeding to the discussion of the map, we review the standard results for volume 



preserving flow in <:5.1 The time \fE flow of this vector field will be shown to approximate the 



map in ^5.2 



5.1 Saddle-Center-Hopf Bifurcation 

For a system of differential equations, the saddle-center-Hopf or Gavrilov-Guckenheimer bifurca- 
tion is the codimension-two bifurcation that occurs when an equilibrium has simultaneously one 
zero and one pair of imaginary eigenvalues. The three-dimensional, center manifold reduction 
to normal form is discussed in depth in |DI96| IDIKOH IGH02| IKuz04j . The key simplification 
is that the "formal" normal form has cylindrical symmetry to all orders in the power series 
expansion: it exhibits the symmetry of the linearized system. This symmetry is typically only 
formal and is broken by terms "beyond-all-orders" . For a system of divergence-free differential 
equations this bifurcation is co dimension-one, and its unfolding is considerably simpler [ BroSl] . 
In the neighborhood of an equilibrium point with eigenvalues (27ria;, — 27riu;, 0), the normal 
form can be most easily obtained in the complex coordinates, (u, v = u,z), that diagonalize the 
linearization. The flow of the linear system then commutes with the rotation u >-^ uexp(2TTi(p). 
Consequently, in symplectic, cylindrical coordinates (r, 6, z) with 

,27Tie „ /^„-27Tie 



u = V2re'"*^ , V = y/2re-'^''' , (9) 

the linearized vector field has the form 

V = uode . 

The normal form has the same symmetry and since the rotation axis r = is a fixed set of the 
symmetry, it is invariant under the dynamics. The normal form vector field therefore has the 
form 

V = rF{r, z)dr + n{r, z)de + Z{r, z)d^ , 

which is divergence free when dzZ -\- dr{rF) = OJj The divergence-free condition implies that 
the two-dimensional projection of vector field onto (r, z) is Hamiltonian with 

H{r, z) = rG{r, z) , 

such that F = —Gz and Z = G + rGr{\ Thus the normal form is a skew-product system of the 
form 

dH 

oz 

e = n{r,z), (10) 

dH 

^= ~^ ■ 
or 



^ The divergence in the new coordinates looks Euchdean because the transformation to symplectic cylindrical 
coordinates is volume preserving . 

^ Indeed, any three-dimensional, volume- preserving flow with a Lie symmetry has an invariant that gives rise 
to a Hamiltonian structure on the projection of manifold by the group orbits |HM98) . Specifically suppose V has 
the Lie symmetry W , i.e., \V, W] = 0. Then if both V and W preserve the volume form Q,, we have dH — —iyiw^ 
and the symplectic form is cj = —iw^- 



Now consider the dynamics of (10) near the origin. The first few terms in a power series 
expansion of the Hamiltonian are 

H{r, z) = r(Ao,o + ^1,0^ + ^o,i-z + ^0,2^^ + ••.)• 

When Ao,2 7^ 0, the imphcit function theorem imphes that the coefficient ^0,1 can be ehminated 
by an affine shift in z. In the new coordinates, we replace j4o,o by —b\ this represents the 
unfolding parameter. The shape of the contours of H near the origin for small b are determined 
by the privileged scaling r = 0{S) and z = 0{yS). In this case, as (5 ^ -ff is equivalent to 

F(r, z) = r(-(5 + /3r + az^) + 0{b^!'^) . (11) 

When a(ib 7^ 0, the shape of the contours of H depends only upon the two signs 

s\ = sgn((5a) , S2 = sgn(a/3) . (12) 

Indeed, the transformation 

1 

lit ^ h = \/\6a\ sgna , 



r, C 



leads to the scaled Hamiltonian 



H{p,C) 



z , r 



S2 



-•51/5+ yP' 



+ pC + 0(\/5) 



(13) 



(14) 



which has no continuous parameters. The implication is that there are two types of saddle- 
center-Hopf bifurcation, see Fig. [6j When S2 = 1, three equilibria are created as ba changes 
from negative to positive. Two, at (p, C) = (0, ±1), are saddle-foci of the three-dimensional flow 
when Q is nonzero. The third equilibrium is a center at (1,0); this corresponds to an elliptic 
invariant circle in M^. 

By contrast, when si = —1, the bifurcation is subcritical: the saddle-foci exist when s\ < 0; 
they annihilate at 6 = 0; and a hyperbolic invariant circle is created for si > 0. 

The separatrix of stable and unstable manifolds of the saddles is the contour H = 0; it is 
the C-axis together with the parabola p = 2s2 (si — C'^). When si = S2 = 1 this separatrix 
encloses a region of closed trajectories that become, for the three-dimensional flow, a family 
of two-dimensional invariant tori that surround the invariant circle. In the original, unsealed 
variables, the volume of the vortex-bubble as delimited by the invariant manifolds of the fixed 
points on the symmetry axis can be easily computed to be 



Vh 



dr Ad9 A dz 



H<o 




(15) 



However if si = S2 = — 1, then the stable and unstable manifolds of the saddles are unbounded 
and there is no heteroclinic connection apart from the ("-axis. 

On any invariant torus, the dynamics is conjugate to a rigid rotation. Since the flow of H 
is integrable, there exist angle-action coordinates (</>, /) that are valid inside the separatrix. In 
these coordinates the normal form vector field would become 



10 




Figure 6: Contours of the Hamiltonian ( 14 1 for the four choices of the signs si = sgnSa and S2 = sgnaf3. The 



upper two panes illustrate the supercritical case, S2 = 1, and the lower two the subcritical case S2 ~ —1. Stable 
and unstable manifolds of the saddle equilibria are shown in blue and red, respectively. 



where z^(/) = dH/dl is the Hamiltonian frequency. Thus the dynamics of the cyhndrical angle 
can be trivially solved to obtain 



e{t) = e{o) + 



1 



m 



n{(t>,i)d4>. 



Consequently, the winding number on the torus is 

As we will see below, the flow is a good model of the map for 6 <^ 1. Consequently, the 
structure that we observed in Fig. [5] corresponds to the creation of a bubble of bounded orbits 
in a supercritical saddle-center-Hopf bifurcation. 

5.2 Saddle-Center-Neimark-Sacker Bifurcation 

Consider a three-dimensional volume-preserving map with a fixed point whose multipliers are 



2TTiul 



Ai = A2 = A = e^-*" , As = 1 , 



(16) 



11 



where, without loss of generahty, a; G (0, j)- As for the flow case, it is convenient to use 
complex coordinates to diagonalize the linearization: let C = (n, v, z) where u is the complex 
eigen-coordinate for Ai, t> = u for A2, and z is the real eigen-coordinate for A3. 

The dynamics in the neighborhood of this fixed point can be studied by a standard normal 
form analysis |Mur03| IBC93| . We summarize the aspects that are new for the volume-preserving 
case in Appendix [Dj 



As for the flow case, the normal form, (43), commutes to all orders with the phase shift 
u -^ ue^^. Consequently, the dynamics of |u| and of z are independent of those of the argument 
of u. To make this explicit, we use the symplectic, cylindrical coordinates (|9| so that if n = x+iy, 
the volume form obeys 



Vl = dx /\ dy /\ dz = -du A du A dz = dr A d9 A dz . 



(17) 



To cubic degree, the volume-preserving normal form obtained in Appendix |D| (43) is 

r' = r (1 - 2az - (7 -h 2a/3)r + (4a^ - 3k)z^) , 
9' = 9 + n{r,z) , 

z' = -6 + z + az'^ + 2(3r + 2jrz + kz^ , 



(18) 



where 



27ri7(r, z) = 2-kuj + AiZ + 2Cir + {aAi + Bi)z^ + (2aCj - 2AiCr)rz . 



Here 4(7,- = —7 — 2a/3, but the remaining eight parameters (q, /?, 7, 5, k, Ai, Bi, Ci) are indepen- 
dent. We summarize this result as the theorem: 

Theorem 1. Suppose fp is a family of C^ volume-preserving maps ofM.^ and that /o has a fixed 
point X* for which 

• DfQ{x*) has eigenvalues (e^'^*'^, 



g-2™ ;l1 



the vector ^ 



g — ip=o is not in the range of Dfo^x* 



Then for p and x — x* sufficiently small, there is a coordinate change so that fp reduces to (18) 
through cubic order. 



The {r,z) components of the map (18) are independent of 9] consequently, the (r, z) map 
is area preserving to 0(3). A typical phase portrait for this two-dimensional map is shown in 
Fig.[7l 



As (5 — > 0, the (r, z) projection of (18) is approximately equivalent to the time h map of the 



flow of the Hamiltonian (14). To see this, introduce the scaled coordinates (13) to the map to 
obtain 

p' = p- 2hpC + 0{h^) , 

(' = c + h{-si + e + S2P) + Oik'') . 



This, to 0{h), is the flow of (14); consequently, the analysis of ^5.1 implies that the phase 
portrait of the map will asymptotically approach the pictures shown in Fig. [6] as (5 — > 0. 
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Figure 7: Phase portrait of the map ( 18 1 for 5 = 0.2, a = /? = 1, and Cr = 7 = k = in the (r, z) plane. The 
unstable manifold of the lower saddle is shown in red and the stable manifold of the upper saddle in blue. Since 
Cr has been set to zero, the invariant circle is artificially a,t z = 0. 



To obtain the portrait in more detail, we study the map directly. Since we have expanded 



in a power series, the valid fixed points of (18) emerge from the origin when 6 = 0. Assuming 



that a/3 j^ 0, there are three such fixed points: 

(0,z±)= (0,± 
{re, Zc) = 




(19) 



where Cr is given by (44) 



The stability of the fixed points can be classified by computing the "residue" 

1 



R 



(2 - tr Df) = a(3r - a'z' + 0(3) 



Recall that fixed points are hyperbolic saddles when i? < 0, elliptic when < i? < 1 and 
hyperbolic reflection-saddles when K> \. 

The fixed points on the z-axis exist when (5a > and have residue 

R^ = -5a + 0{6^''^) . 

Thus these points are always hyperbolic saddles when they exist. The z-axis corresponds to 
the one-dimensional stable/unstable manifolds of these two fixed points with the corresponding 
multiplier A-t = 1 -|- 2az± + 0{5^''^). Thus when a > the z-axis is the unstable manifold of the 
upper fixed point and the stable manifold of the lower fixed point. For the three-dimensional 
map, the other two multipliers form a complex conjugate pair, 

A± = (1 - 2az±)e2'^*^± + 0{6^/^) . 

provided that uj± = 0(0, z-t) ^ 0. The graph of two-dimensional invariant manifolds for these 
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fixed points lias the for 



»'('■' = ^-- 2^1^ -d^'-'+°(^'' 



Note that the manifold of the upper (lower) fixed point has negative (positive) slope when 
a/3 > 0; these manifolds are thus inclined towards each other and, as we observed in Fig. [TJ 
they generally intersect, enclosing, roughly speaking, a topological sphere. When af3 < the 
manifolds are inclined away from each other and, as 5 — > 0, there is no local intersection of the 
manifolds. 



The third fixed point corresponds to an invariant circle C of the three-dimensional map (18). 
Since r > 0, the circle exists only when 50 > 0. The dynamics restricted to the circle is a rigid 
rotation with rotation number 

UJL = ^{Vc, Zc)=UJ+ ^r^ {aCi + CrAi) 5 + 0(6^) . (20) 

znap 

We call this the longitudinal rotation number of C. 

The invariant circle corresponds to a fixed point in the (r, z) plane with residue 

Consequently if a/3 < 0, the invariant circle is transversely hyperbolic. Alternatively, if a/3 > 
the invariant circle is elliptic when a6 < 2. Nearby circles in the (r, z) plane have a rotation 
number defined through R = sin^(7ra;T), 

-T = ^ + 0(5^^) . (21) 

This is the transverse rotation number of the invariant circle C. 

5.3 Quadratic Map near the Saddle-Center-Hopf Line 

When e is small and —4 < // < 0, the normal form (IT]) can be transformed to the saddle- 



center-Hopf form, (18). As this transformation is somewhat tedious, we use computer algebra 
to perform the manipulations. 

The fixed point at the origin for e = and — 4 < ^u < has the linearization 




To apply the results of ^5.2 this matrix is transformed to the diagonal form M, (38), using 



1 


1 


1 


M 


M 




A 


1 





A-1 


A-1 


1 


1 






^ To make sense of this expansion, we should use the scaled variables (|13|). This scales the bubble to fixed size 



as J — > 0. In this case each of the terms in the series begins with O^S'^) terms 
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Applying this transformation to M gives the new map T^^ o f o T. Note that the Jacobian 
of this transformation is det T = i cot ttuj, so the volume in the transformed coordinate system 
must be scaled by this factor. The transformation ip to the new coordinates is constructed by 
Lie series and is obtained first to 0{e^) through second order in ^ and then to first order in e. 



As discussed in Appendix [Dl this gives a map of the form (47). A final affine transformation on 



z restores the magnitude of the multipliers to 1, giving a map of the form (18) with parameters 



a ^ 2a + {h-2c)ii 2 a - 6 3 + // K, ,: e 

a = , P = 3 , 7= j-r, — ^ /^' ^ = 



Ai 
Bi 
Ci 



a — b 



{4{n + 2)(^ + 3)a^ + /i(4(^ + 4)c - (5/i + 16)b)a + b^fijfi + 2)) 
2M-Mm + 4))3/2 ' 

-2{ii + l)a2 + //(2c(5// + 13) - 6(3^ + 5))a - /i (4/x(2/i + 5)c2 - 26//(3/x + 7)c + b'^ifi + 1)^) 



47ra^Y^-/i(4 + ^u) 






(22) 



where iOo{fJ-) is given in n. Note that these coefficients have singularities when /i G {—4, —3, 0} 
corresponding to the resonant frequencies a; e {^, ^,0}, respectivelyrj 

Thus the dynamics of (II]) near the saddle-center-Hopf line are approximated by the normal 
form (18) away from the low-order resonances. The signs (12) that determine the character of 



the bifurcation now become 



si = sgn((5a) = sgn(eo) , 

S2 = sgn(a/3) = sgn(a) sgn(2a + (6 — 2c) /x)) . 



Recall from ^5.1 that the bifurcation is super(sub)-critical for S2 = 1 (= ~1)) see Fig. km For the 



parameter values {a,b,c) = (1, |, ^) that we have primarily studied, S2 = 1 for all ^ < 0, thus the 
bifurcation is always supercritical. Consequently as e increases through zero, an elliptic invariant 
circle and a pair of saddle-focus fixed points whose two-dimensional manifolds intersect are 
created. This is confirmed by numerical computations; for example in Fig. [5] where e = 0.01, the 
map is nearly integrable and the ball enclosed by the two-dimensional manifolds is predominantly 
filled with invariant tori. As e increases the intersection angle between the two-dimensional stable 
and unstable manifolds of the saddle-foci grows, many of the tori are destroyed, and orbits near 
these manifolds become unbounded, see Fig. |8j Along the way there are many torus bifurcations 
that create new families helically wound around the original invariant circle. We postpone a 
discussion of these until f|7j 

Meanwhile, these results permit comparison with Fig. [4] near e = 0. Taking into account the 
scaling of the volume due to the transformation T and that due to the complex transformation 



(17), the fraction of bounded orbits in the box B, mn is 



Vh vra 4: + fj, 

'vol(e)detr " 96/3 V^^ 



^rei = 2^;;;;T7RY:ii7^ = ^^\/3:;^' (^3) 



* The value {b — 2c)^ + 2a = is also special since at this point both 13 and 7 vanish. This corresponds to the 
boundary between supercritical and subcritical saddle-center-Hopf bifurcations. To study the dynamics at this 
point would require computing higher order terms. 
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where Vh is the volume (15). This formula gives the dashed curve shown in Fig. El Since 
P = 0{fi~^), (23) implies that the bounded fraction is 0(^^'^) when /i — > 0~. The actual results 
in Fig. H] deviate from this result when |/i| < 0.05 due to the a; = resonance. Moreover, the 
resonances near fj, = —2 and —3 are also not captured by this lowest-order formula. However, 
(23) agrees well with the numerical results near the period-doubling point /i = —4. 





Figure 8: Orbits of the map ^ ior n = -2.4, a = 1, 6 = c = 5, and e = 0.1 (left), 0.2 (center), and 0.32 (right). 
In the last case there appear to be no bounded orbits, however there is probably a pair of period 7 saddles. The 
orbits that are shown lie near the stable and unstable manifolds of these orbits. The red cubes are centered on 
the origin and have sides with length ^\/e. 



From (22) we can compute the position of the fixed points and invariant circle (19). These 



expressions are rather complicated and not especially useful. However, the rotation numbers of 
the invariant circle will be used extensively in the following sections. For the values (a, b, c) = 

1 l^ 

2' 2/ 



(1, i, i) these become 



UJL = tJo(/W) + 



U>T 



5/i3 + 19/i2 - 28/i - 64 
167Tfl{fI + 3)(/i - 4)(-^(// + 4))3/2 



gifj-) 



e + Oie^ 



vr/z V 2 
58/i'^ + 648/i'' - 377^5 - 23526/ - 59027//^ + 139032/^^ + 689616// + 695808 



(24) 



96/x3(/i-4)2(/i + 3)2(/x + 4)2 



As a first step in understanding the invariant tori, we will study the evolution of the invariant 
circle. As we will see, resonances between the rotation numbers (loljUJt) will be key to this 
discussion. 

6 Invariant Circles 



In a volume-preserving flow, the supercritical, saddle-center-Hopf bifurcation creates a pair of 
saddle-focus fixed points, one of type-(2, 1) and one of type-(l,2) separated by 0{\/^), recall 



^5.1 The two-dimensional unstable and stable manifolds of these fixed points coincide in the 
formal normal form and bound a ball that is foliated by a family of two-dimensional invariant 
tori. These tori enclose an invariant circle that has a diameter of size 0{y6). 



As we showed in ^5.2, the analogous bifurcation for a volume-preserving map has a normal 



form (18) that is approximated by this flow when 6 <^ 1. However, the cylindrical symmetry of 
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both normal forms is only formal: even when 5 is very small, the invariant circles presumably 
exist only for a Cantor set of parameter values when their rotation numbers satisfy appropriate 
Diophantine conditions fCS90a| IXia96j . 

Our goal in this section is to study numerically the invariant circles for the quadratic map 
M that are created in the saddle-center-Neimark Sacker bifurcation. In order to do this, we 
must first develop an algorithm to accurately compute elliptic invariant circles. 

6.1 Ellipsoid Algorithm 

A number of algorithms have been proposed for finding invariant circles and tori; some of the 
history is discussed in |SOV05| IHDLL06b] . Many of the techniques assume that the dynamics on 
the invariant circle is conjugate to a rigid rotation with a given rotation number uj^. In this case, 
the conjugacy can be expanded in a Fourier series and a Newton method employed to compute 
the Fourier coefficients. Castello and Jorba have used this idea for Hamiltonian systems that 
have a Cantor set of circles whose frequencies vary |CJ00j . If the gaps in the Cantor set are 
"small" they can fix a rotation number and search for the corresponding circle. Haro and de la 
Llave |HDLL06al IHDLLOGbl IHDLL07] also use spectral methods to find tori in quasiperiodically 
forced maps. The important simplification for the forced case is that the rotation vector and the 
conjugacy to rigid rotation are automatically known. A very different method was developed by 
Edoh and Lorenz |EL03] to find nonsmooth, but attracting invariant circles; here the asymptotic 
stability of the invariant circles is key. A related method is to define a variational method, using 
the distance between the circle and its image, and obtain the circle by following the gradient 
flow jLCCOGj . 

An iteration-based method was developed by Simo |Sim98| . Letting T = — be the "irrational 

period" , Simo defined the Tth iterate of the map by interpolation between the points / L -I and 
yr^l J5J An invariant circle corresponds, roughly speaking, to a fixed point of /^. An advantage 
of this method is that it does not depend upon the existence of a conjugacy. It can also be used 
to determine ujl by fixing a section and looking for approximate returns to define T |CJ00j . 

Our method is also iteration-based, but in contrast to Simo's idea, uses a slice to localize 
iterates and a number of points to locate the invariant circle. Suppose C is an elliptic invariant 
circle with a local cross section S; for ([T|, the section 

S = {(x,y,0):y>0} 

works in most cases. Let Sa = {\z\ < A,y > 0} be a thin slice enclosing S. Assume that 
there is a neighborhood of C in which the local dynamics predominantly lies on two-dimensional 
invariant tori that enclose C. Thus the intersection of a typical orbit near C with Sa is a slightly 
curved, elliptic cylinder, see Fig. [9j 

Begin with a point (^o £ 5] in the neighborhood of C and find a set of A^ returns of its orbit 
to Sa: 

C = /*HCo)eSA, 0<h<t2<...<tN , i = i,...,N. 

By assumption, these points lie on an invariant torus; the goal is to find the "axis" of this torus. 
Instead of fitting the points to an elliptic cylinder, it is easier to fit them to a general ellipsoid, 
defined as level set of a quadratic form 

e{x, y, z) = Ax^ + By^ + Cz"^ + Dxy + Exz + Fyz + Gx + Hy + Jz . 



Note that these points may be far apart and so the interpolation may not be accurate. 
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Figure 9: Finding an invariant circle by fitting a section of a nearby torus to an ellipsoid. 

Since the returning points are, presumably, all close to the circle and since they lie on the 
ellipsoid by hypothesis, consider the A^ differences 5C = C ~ Co- In this coordinate system 
the ellipsoid goes through the origin, so e(5C*) = 0. Since all the 5^^ approach zero when the 
method converges we normalize by / = -^ X^j=i l<^C*l a-iid introduce 6^ = ^C /^- This removes 
the singular behavior from the algorithm that originates from the fact that the fitted ellipsoid 
becomes thinner, and it also helps avoid round-off errors. 

Since the equations e((5^*) = are homogeneous in the coefficients, an additional equation is 
needed to fix the solution; we chose to normalize the sum of the squared coefficients of e, giving 
A^ + 1 linear equations in the coefficients: 

e{6e)=0, i = l,...,N, 
A + B + C = 1 . 

We must choose at least A^ = 8 to fix the 9 coefficients of e. It seems most convenient to use 
exactly 8 returns, since the number of iterates to find a return can be quite large, especially 
when the slice thickness. A, is smalljj 

Once the ellipsoid is determined, the center of the ellipse in the section S is used as the next 
guess for a point on C; for the section S = {z = 0}, this gives the iterative step 

Co ^ Co + ^ ^aJ_b2 (^H - 2CG, BG - 2AH, 0) . 

The whole process is then repeated, finding a new set of returns to the slice Sa and fitting 
a new ellipsoid. The error is estimated as the change in distance between the new point and 
the previous one. The iteration stops either when the error decreases below an error tolerance 
(we typically chose 10"^*^), or after a fixed number of iterations (we chose 10). If the given 
tolerance is not obtained for a fixed A, we decrease A by a factor of 100 and then restart the 



* This could be a problem if the points do not lie in a general position; to get around this one could use more 
points and solve the system using least squares. 
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algorithm with the best previous value as the initial point Co (the results typically converged 
when A = IQ-^). 

Indeed, it is not hard to 



This process appears to be quadratically convergent, see Fig. 10 



show that this is the case for a two-dimensional version of this algorithm for finding an elliptic 
fixed point: since the linearized map has invariant ellipses the linearization of our iteration at a 
solution is superstable. 

We used our algorithm to compute invariant circles for the normal form (II]) for a grid of 
parameters near the saddle-center-Hopf line with the standard values a = l,6 = c=2- Typically 
we started at a fixed value of /x with e = 0.01. If we are able to find the invariant circle, e is 
incremented using the previous solution as an initial guess for the new parameters. 
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Figure 10: Error on a log(log) scale for successive iterates of the invariant circle position for fi = —2.4, with 



10 and the values of e indicated. The line represents an error of e 



which would mean quadratic 



convergence. With this value of A, the minimal error seems to be about 10 ^^ , and is achieved in 4 — 6 iterations. 

An example of the output of this algorithm with ^ = —2.4 is given in Table [T] For this 
value of 11 the algorithm finds an invariant circle when e < 0.3129. Beyond this value it no 
longer converges; moreover, when e > 0.326 the computation of §3 (with a 500^ resolution and 
10^ iterates) also gives no bounded orbits. Consequently we believe that the circle is either 
destroyed or becomes unstable in this range of e. Several of the computed invariant circles are 
overlaid in Fig. 11 This figure also shows a chaotic orbit for e = 0.315 near where the invariant 



circle used to be. 

The ellipsoid algorithm can fail in two ways. One the one hand if the invariant circle has a 
very small stable neighborhood, the initial guess Co may lie on an unbounded trajectory even 
when there is an elliptic circle. This typicially happens as the parameters approach a bifurcation 
where C becomes hyperbolic. To approach the stability border, we can simply take smaller steps 
in e. This is apparently what happens for e > 0.31 when /i = —2.4. A second failure mode 
occurs when the orbits remain bounded, but the given error tolerance cannot be achieved; this 
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e 


X 


y 


error 


UJL 


Wt 


0.01 


-0.000077100956 


0.110334966059 


1.6e-ll 


0.282171317669 


0.00939482 


0.02 


-0.000218667148 


0.156587932592 


2.7e-ll 


0.282294227518 


0.01331187 


0.03 


-0.000401253015 


0.192348589289 


7.8e-12 


0.282415802529 


0.01633514 


0.04 


-0.000606226180 


0.222722556427 


1.5e-ll 


0.282536006132 


0.01889745 


0.05 


-0.000263469053 


0.250886913814 


l.le-11 


0.282654815793 


0.02117252 


0.06 


-0.001281607227 


0.273786477131 


1.4e-13 


0.282772154554 


0.02324300 


0.07 


-0.001578623424 


0.296421783954 


8.4e-ll 


0.282888017858 


0.02516009 


0.08 


-0.001928385257 


0.317549729412 


9.6e-12 


0.283002357933 


0.02695743 


0.09 


-0.002313763896 


0.337472296148 


3.3e-12 


0.283115135201 


0.02865947 


0.10 


-0.002732253247 


0.356389819268 


7.2e-13 


0.283226310896 


0.03028362 


0.11 


-0.003191438702 


0.374443365374 


l.le-11 


0.283335847312 


0.03184166 


0.12 


-0.003632379081 


0.391797224416 


8.9e-ll 


0.283443708609 


0.03309808 


0.13 


-0.004150372053 


0.408459677586 


3.4e-12 


0.283549859693 


0.03480205 


0.14 


-0.004689339128 


0.424551322589 


2.5e-12 


0.283654271356 


0.03622014 


0.15 


-0.005256682085 


0.440125047284 


2.9e-ll 


0.283756917316 


0.03760506 


0.16 


-0.005618175619 


0.455286683823 


5.7e-ll 


0.283857779176 


0.03896231 


0.17 


-0.006487124526 


0.469903756471 


5.7e-12 


0.283956842148 


0.04029660 


0.18 


-0.007146262564 


0.484186136060 


4.7e-12 


0.284054109448 


0.04161036 


0.19 


-0.007838687257 


0.498105595397 


6.6e-12 


0.284149593962 


0.04291542 


0.20 


-0.008559583750 


0.511685650160 


4.9e-ll 


0.284243328335 


0.04420893 


0.21 


-0.009315347590 


0.524951884960 


6.6e-ll 


0.284335369737 


0.04549702 


0.22 


-0.010104271994 


0.537923761008 


2.9e-12 


0.284425807465 


0.04678551 


0.23 


-0.010926038044 


0.550616855400 


5.2e-ll 


0.284514773462 


0.04807963 


0.24 


-0.011781835115 


0.563047279372 


1.4e-ll 


0.284602457095 


0.04938499 


0.25 


-0.012671957974 


0.575230413835 


3.6e-12 


0.284689126532 


0.05070952 


0.26 


-0.013597339065 


0.587176221467 


9.2e-ll 


0.284775161033 


0.05206063 


0.27 


-0.014558693626 


0.598897558338 


2.4e-ll 


0.284861102396 


0.05345099 


0.28 


-0.015555864549 


0.610403755970 


2.8e-13 


0.284947742217 


0.05489500 


0.29 


-0.016590212710 


0.621704492198 


5.2e-12 


0.285036283041 


0.05642285 


0.30 


-0.017662158192 


0.632808106481 


3.0e-ll 


0.285128671076 


0.05807352 


0.31 


-0.018772488310 


0.643722521310 


8.1e-12 


0.285228402490 


0.06007135 



Table 1: Position (j;,y,0) on S and rotation numbers of the elliptic invariant circle of (IT]) for 
/i = —2.4, and (a, b, c) = (1, ^j 5) ^^ function of e. The error is the distance on S between the 
last and the penultimate iterates. We attempted to find a circle with an error < 10~^^ using a 
slice thickness A < 10~®. 
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Figure 11: Several of the invariant circles for /i = —2.4 and e values as shown. Each circle is displayed with 
1000 iterates. For the outermost set of points, at e = 0.315, there is apparently no invariant circle. 

indicates that there is no elHptic invariant circle. For example, when ^ = —2.4, the algorithm 
does not converge for 0.3018 < e < 0.3044; in this range, the circle appears to be unstable and 
to have undergone a doubling bifurcation, see ^ 

In our computations for other values of /i, we observe that the invariant circle often undergoes 
a complex sequence of bifurcations as e changes, sometimes loosing stability or simply vanishing 
for one e and then regaining stability or reforming for a larger value of e; more details are 
discussed in ^ below. In every case there appears to be a maximal value of £ beyond which 
the circle ceases to exist or never regain stability. To make more sense of these bifurcations, we 
compute the longitudinal and transverse rotation numbers. 

6.2 Longitudinal Rotation Number 

If there is an invariant circle, then the restriction of the map to the circle, f\c, is a homeo- 
morphism. Recall that the rotation number of a circle homeomorphism always exists and is 
independent of the initial point. This is the longitudinal rotation number of the circle, lvl- 

The longitudinal number can be easily computed under the assumption that the invariant 
circle projects to the plane j; = as a Jordan curve that encloses the the originjj and that the 
angle advance between iterates "avoids an angle" |DSMOOJ . The avoided angle is chosen as the 
location for the branch cut of the arctan function. Typically the avoided angle can be chosen as 
the angle diametrically opposite to the average rotation angle. Then we simply sum the polar 



"^ If this is not the case, we could also use the self-rotation number |DSM00| . 
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angular increments between images, defining 

1 ^ 
G){N) = --y^atan2{ytyt-i + ztZt-i,ytZt^i - ztyt-i) 



27r 
t=i 

where atan2(y, x) = arg(x + iy). An approximate rotation number is then simply Q{N)/N + 

O(iV-i). 

However as was first suggested by Henon |EV01] . a more accurate value for wl is easily 
obtained by carefully choosing N . To do this, compute the continued fraction expansion of 
ujl by constructing the sequence, — , of its continued fraction convergents. Each convergent is 
defined to be a rational number closer to loi than any others with denominators q < Qj- Since 
each orbit on C must be ordered as a rigid rotation with rotation number ujl, we can obtain the 
convergents by finding the closest returns to the initial point. Let qj be the time of the next 
closest return to the initial point, dj = \Cq^ — Col < dj-i, and pj = [@{qj) + gj be the nearest 
integral number of rotations. To start the process we arbitrarily chose dj^ = 0.01; this means 
that the first computed convergent will not be the leading convergent of the continued fraction 
for ujl; however, it is easy to find the earlier convergents from the continued fraction expansion 
for pi/qi- For example with (e, ^) = (0.2, —2.4), the computed convergents are 

qj = {971, 3205, 7381, 150825, 158206, 309031, 776268} , 
Pj = {276,911,2098,42871,44969,87840,220649} , 

which implies that the continued fraction for ujl is 

220649 

UJL ~ = [0, 3, 1, 1, 13, 3, 3, 3, 2, 20, 1, 1, 2] . 

776268 L 7 , , , , , , , , ) 1 > J 

Note that this rotation number is relatively close to the rational i = [0,3, 1, 1] — see §7 Recall 
that the rotation number is close to its convergents: 



Pj 

UJL- — 

Qj 



1 
< — 



Indeed, this is the error observed in practice: Fig. [12] shows the difference between the value 
at the jth step and the final value. As discussed in [EVOlJ . a value correct to 0{q~ ) can be 
obtained by adding to ujl the average angular deviation from exact q'j-periodicity over the next 
qj iterations; however, the estimate above is sufficient for our purposes. 

The values in Table [l] were computed using the sequence of closest returns for t < 10^, 
giving a result accurate to 0(10"^^). We can also use this computation to check that the orbit 
is properly ordered on the circle. Indeed, the sequence qj must obey the recursion 

Qj+i = (^jQj + Qj-i ' 

where q-i = 0, and qo = 1, and Uj € N are the continued fraction elements. Thus a necessary 
condition for a valid sequence of convergents is that 

{qj+i - qj~i) mod qj =0 . 
Usually this criterion only fails when the error bound in the circle algorithm is also not achieved. 



The computations of u>l are compared with the normal form results (24) in Fig. 13 The 



dominant behavior is the zeroth order rotation number uJo{fi) in ([7|. For the figure we subtracted 
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Figure 12: Error in tlie computation of tlie rotation numbers of an invariant circle at /x = —2.4, for two different 
e values as shown. The open points correspond to oj^ and the solid to lot- The error in the longitudinal rotation 
number is 0(9^^), but that in the transverse rotation number is only 0(g~^). 

this value from the computed results and then compare with the theoretical 0{e) term. The 
curves in the figure show the results for fixed values of e as a function of [i. The agreement 
between the numerical results and the theory is nearly perfect for e < 0.15 away from the 
resonances where the normal form is not valid. The computations show that an ellipitical 
invariant circle does not even exist near the main resonances /i = —4, —3, —2 and 0, see ^ The 
numerical results indicate an additional singularity in ojl near /i = —2 that is not present in the 
normal form to 0(e). To find this singularity in the normal form we would have to keep quartic 
terms. 



6.3 Transverse Rotation Number 

Assume that / does have an invariant circle C : {C(^l) '■ Gl £ §^} that is C^ and on which the 
dynamics is conjugate to rigid rotation with irrational rotation number ujl, 



The linearization of / then gives rise to a quasiperiodic skew-product on 

0' = eL+ UJL mod 1 , 



xS^ 



(25) 



where A{9i,) = Df{C,{6L)) is a periodic matrix. 

The transverse rotation number of C is the average rotation rate of a transverse vector ^ 
"around" C, if that average exists. Indeed, Herman proved that the "fiberwise" rotation number 



for a quasiperiodically forced circle map always exists |Her83] . The map (25) reduces to this 
case if the ^ dynamics are projected onto a family of circles transverse to C. To do this, define 



the transverse angle ip relative to a ribbon TZ attached to the invariant circle, see Fig. 14 The 
transverse dynamics then induce a circle map ip 1— > g{ip, Ol)- Consequently, under the assumption 
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Figure 13: Comparison of the computed Aojl = ^l ~ ^o{l^) (dots) with the normal form (24 1 (curves) as a 
function of /i for (a, 6, c) = (1, |, |) and the values of e indicated. 



■*-.., 




Figure 14: Computing the transverse rotation number for an invariant circle C relative to a ribbon 7?.. 



that C exists and its dynamics are conjugate to a rigid rotation, the rotation number of g exists 
and is independent of the initial 9 and <p. 

Let t(^L) be the unit tangent vector to C at C(^-l); it can be approximated using the closest 
approaches from the computation of lol. Since the invariant circle of M appears to be almost 
always everywhere transverse to the x- axis rl define the ribbon by the vector r{9L) = ^(^l) x ^i- 
Beginning with an arbitrary vector vq attached to the point ("(0) G C, we iterate to obtain the 
sequence Vj+i = A{ju!l)vj. These vectors are then projected onto a plane orthogonal to the 



An exception is shown in Fig. 
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for {e,fi) = (0.068,-1.9), where the invariant circle develops "curlicues". 
Here the algorithm computes ujt incorrectly because the projection of the circle on the a; = plane is not 
one-to-one. 
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local tangent vector t{jujL)] for convenience, we also rescale, defining 

p = \r\ [v — {t ■ v)t) . 

The projected vector is effectively two-dimensional with components p = {p" , p ) defined relative 
to the ribbon direction 

ptt = p . f = ei • t X u , 

p = i ■ p X f = ei ■ {v — {i ■ v)i) . 

Thus the angle oi p is ip = arctan^p-^/p"). As before, we compute the change in rotation angle 
at the jth iterate with the two- argument arc-tangent: 

Aipj = atan2(p]_-^p^! + pf-iPf,Pj-iPf - pf-iPj) ■ 
The transverse rotation number is the average change along the orbit 

i=i 

It is not clear how to optimize the error in this computation as we did for the longitudinal rotation 
number. We simply use the time of closest approach that we computed for the longitudinal 



rotation number; however, as can be seen in Fig. 12 the accuracy for this computation is only 



We compare the computations of lot to the normal form results (24) in Fig. 15, Again, the 
agreement between the two results is extremely good for small e and away from the resonant 
values. The C(e^'^) normal form result is finite at the resonances at /i = —4, —3, and —2, but, 
as before, the numerical computations indicate that the actual rotation number diverges there. 
Such divergences are found in the normal form at 0(e^'^); however, since these corrections to 
the rotation number are normally very small, we do not show in them in the figure. 

6.4 Frequency Map 

If /i and {a,b,c) are held fixed, the rotation numbers (culjUJt) of the circle vary along a curve 



as e changes; an example is shown in Fig. 16 for fj, = —2.4. For small e this curve lies close 



to the parabola (24) defined by the normal form results. Though we expect the curve to be 
defined only for a Cantor set of parameter values, it appears to be continuous for most values 
of e; there are, however, several small intervals visible in which the ellipsoid algorithm does not 
converge. The algorithm can find an invariant circle up to e = 0.312, at this point the circle is 
apparently destroyed by a resonant bifurcation, see ^ The invariant circles shown in Fig. 16 



are also discussed in more detail in ^ 

For some values of n, the frequency map exhibits singularities as e varies; two examples are 



shown in Fig. 17 Note that the normal form curve (brown) still fits the results for small e. The 
singularities are in col; those shown occur near {e,n) = (0.1,-1.5) and (e,/u) = (0.07,-1.9). 
Near the singularities, the original circle (blue) becomes highly distorted and as e increases it can 
no longer be found. The singularities appear to be associated with a pair of circle-saddle-center 
bifurcations: the old circle is destroyed in a saddle-center bifurcation below the resonance line. 



while a new elliptic invariant circle is born above (red), see ^7.3 Remarkably, the new circle 



persists to a much larger value of e and its rotation numbers approximately follow those of the 

25 




Figure 15: Comparison of the the measured transverse rotation number (dots) with the theory ( 24 1 as a function 
of fi for (a, b, c) = (1, |, |) and the values of e indicated. 



normal form. In each case the old circle eventually appears to lose smoothness and is apparently 
destroyed. 

The image of the frequency map for /i = — 3.9 is shown in Fig. [T8| For this // the longitudinal 



rotation number decreases with e as predicted by (24). The computed curve is fairly smooth, 
but there are bifurcations along the way, e.g., for e = 0.0767, but the resolution is too rough 
to resolve them. The gap between 0.0173 < e < 0.0177 is caused by a period doubling and the 
doubled circle is shown in green. 

A summary of our computations is given by the two-parameter frequency map 



shown in Fig. 19 Here we computed the invariant circle for /x G [—4,0] in steps of 0.1 and for 



e > 0.01. We were able to find invariant circles for some range of e for each /i except for the 
values {0, —2, —3, —4} where there are strong resonances. When —0.6 < /i < 0, an invariant 
circle is stable only for e < 0.01, so these points do not appear in the figure 

There are many cases in which the algorithm fails to find a circle at one value of e, but 
then succeeds for a slightly larger value. These gaps are typically very small, but can be seen 
in several of the fixed /j, curves; in most cases they correspond to resonant bifurcations that 
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Figure 16: Numerically computed frequency map for the invariant circle with ^ = —2.4 and the standard values 
(a, 6, c) = (1, |, I) as e varies from 0.01 to 0.312 (black dots). The curve (brown) represents the normal form of 



( 24 1 . Insets show several phase portraits of the circle and/or nearby orbits. The illustrated resonances in the 
order of increasing e are (7, 1, 2), (3, 4, 1), (3, 3, 1), (10, 3, 3), (46, -2, 13), (7, 0, 2). 

either cause the circle to lose transverse stability or to be destroyed. It is to the study of these 
bifurcations that we turn next. 



7 Resonant Bifurcations of Invariant Circles 

The rotation vector oj = {ujl^ojt) is resonant if there are integers {mi,m2,n), not all zero, such 
that 

m ■ uj = n . (26) 

If Lu does not satisfy such a relation then it is "nonresonant." 



Since (26) is homogeneous, if {m,n) is a solution then so is {lm,ln), for / G Z; consequently, 

(27) 



the set of integers that satisfy (26) form a sublattice of Z^, 

C{u)) = {{m, n) G Z : m ■ uj = n} 



called the resonance module. The dimension of C is the number of independent resonance 
conditions. Using homogeneity, we will assume that n is nonnegative and gcd{mi,m2,n) = 1. 



In w-space, (26) defines a line for each {m,n); a few such lines are shown in Fig. |19[ Using 



the normal form results (24), the resonance curves for the invariant circle can also be obtained in 



parameter space. The fine structure in the bounded volume Fig.[3]is caused by these resonances. 

For example, resonances with \m2\ = 2 are often observed to result in the destruction or loss 

of stability of the invariant circle C, as well as a dramatic drop in the volume of bounded orbits. 



Figure 19 shows that C appears to lose stability at the (1, —2,0) resonance for —1.4 < /i < 
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Figure 17: Numerically computed frequency map for ^ — —1.5 and —1.9 (dots) and the corresponding normal 
form results (brown curves). Also shown are the computed invariant circles for several values of e. The singularities 
in the frequency maps are due to (5,-1,1) and (4,1,1) resonances that result in pairs of circle-saddle-center 
bifurcations, see 37.31 



and at the (3, ±2,1) resonances for —3.3 < n < —2.5 . In each case, a new elliptic, doubled 
circle is created when the original circle crosses the resonance; phase space portraits are shown 
in the top and bottom panels of Fig. 20 Since ujt/^l = 2 ^°^ ^^^ (-'-' ~^' 0) resonance, tori near 



C just below the bifurcation approximate a Mobius strip with a half-twist, and the new circle 
undergoes one-half turn about the original circle for each longitudinal period. For the second 
case, the (3, ±2, 1) resonance, the new circle rotates by 3 turns transversally in 2 longitudinal 
rotations. In both cases the new elliptic invariant circle is a double covering of the original one, 
as is familiar for period doubling of periodic orbits of flows. 

Interestingly, there is also another type of period doubling bifurcation that cannot occur in 
flows. This happens, for example, at the (4,2, 1) resonance; a phase space portrait is shown in 
the middle panels of Fig. [20} Instead of one invariant circle that is a double cover of the original, 
here there are two disjoint invariant circles, each covering the original circle once. Though 
the new invariant circles are geometrically disjoint, they are mapped onto one another by the 
dynamics. The reason for the different behavior is that gcd(?TT,i, 771,2) 7^ 1 in this case. This is 
the simplest example of a general scheme that is easily understood by studying an integrable 
map of the two-torus with resonant rotation vector to. 



Lemma 2. The orbits of the torus map 

x' = X + uj mod 1 
for x,uj (z T^, are either 

1. dense on T'^ if uj is nonresonant (AmiC{uj) = Oj; 



(28) 
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Figure 18: Frequency map for fi = —3.9 as e ranges from 0.01 to 0.173. Insets show invariant circles for three 
values of e. The doubled circle in the gap near e = 0.174 is caused by a (7, —2, 3) resonance. 

2. dense on k circles ifuj satisfies one resonant relation andgcd(mi,m2) = k (dim. C{lo) = 1); 
or 

3. periodic if uj satisfies more than one resonance relation (dim. C{uj) = 2). 

Proof. If UJ is nonresonant, then the result is standard |CFS82j . Suppose now that there is a 
single resonance relation (dini£ = 1), let k = gcd(mi,m2) be the common divisor, and define 

(mi,m2) = k{p,q) , (29) 

where p, q are coprimeP] As is shown in |HW79j there is a solution {p, q) G Z^ to the equation 

pq-qp=l , 

if and only if p, q are coprime. Taking any such solution, define the unimodular transformation 

fp q\ 

\p q 



^ Note that if one frequency is rational, then mim2 — and hence pq — 0. This case can be formally included 
here with k = m\ + m,2 and p or q equal to 1. 
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Figure 19: Frequency map (uLyi^r) for fi G [—3.8, —0.6] with step 0.1 and a grid of e values from 0.01 up to the 
largest e for which the ellipsoid method converged to an invariant circle. Also show are several of the important 
resonance lines defined by (|26||. 



The variables y are simply a new set of coordinates on T^. The map in these new coordinates 
takes the form 



2/2 



n 



y2 + puJi + quj2 . 



(30) 



By assumption k and n are coprime integers, so the orbit of yi is periodic with period k (in 
particular if n = 0, then k = 1). Moreover, since u; satisfies exactly one resonance relation, the 
combination {p,q) • uj is irrational. Thus the orbit of y2 is dense in T^. Hence the combined 
orbit on T^ is dense on k circles. This still holds true in the original coordinate system, but the 
circles will wrap around the torus in both directions. 

Finally, if u; satisfies two independent resonance relations (dim£ = 2), m ■ lo = n and 
m ■ UJ = h, then the components of io must each be rational because then 



m 



, m 



CJ 



and the matrix on the left is nonsingular since the vectors m and m are not parallel. Thus both 
frequencies are rational, and all orbits are periodic. D 

Structures similar to the second case of Lem. |2] are also found in quasiperiodically forced 
circle maps pK06| IJS06J where they have been called {k, (7)-invariant graphs. 



30 



(1. 056,-1.0) 




Figure 20: Circle doubling bifurcations. Shown are some of the orbits in the neighborhood of an ellipitical 
invariant circle (left panels) just before the doubling bifurcations, and orbits in the neighborhood of the doubled 
circle (right panels) just after. The corresponding resonances are (1,-2,0), (4,2,1), and (3,-2,1) from top to 
bottom. The values of (e, fj,) are shown. 

We now embed these dynamics into three dimensions as invariant sets of a family of "inte- 
grable" volume-preserving maps /o on the solid torus, T^ x M"*", 

Ue,r) = {e + u{r,fi),r). (31) 

Here r represents the radius, so that r = corresponds to the toroidal or longitudinal axis, 9l is 
the toroidal angle, and 6t is the poloidal or meridional angle. We assume there is a parameter 
fi to unfold the frequency map (just as for our model (IT]) 
into invariant tori r = ro, and on each such torus /o becomes (28). 

Now suppose that at (ro, /io) 
that is u;(ro, /xq) ■m = n where m satisfies (29). Each orbit on this resonant torus densely covers 
a circle {k 



The phase space of (31) is foliated 



the rotation vector of ^1\ satisfies a single resonance condition. 



1) or set of circles {k > 1) on this torus. These curve(s) correspond to a torus 
knot (torus link) [Ad_a04j ; several examples are shown in Fig. 21^^ A torus knot of type {p,q) 
coprime corresponds to the closed loop {(gt, — pt, ro) : t G M} that has p poloidal wraps for 
each q toroidal circuits. A torus link of type (mi,m2) is the collection of k loops given by 
{(gi, -pt + i/mi, ro) : i = 0, . . . , /c - 1, t G M}. 



^"^ Famous knots and links include the trefoil knot (3,2), the Hopf link (2,2), and Solomon's knot (which is a 
link) (4, 2). Note that Borromean rings are not included in the class of torus links. 
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For example, the (3,4, 1) resonance leads to a single (3,4) torus knot, the (4,4, 1) resonance 
leads to four (1,1) torus knots (each of which is not knotted), which together give the (4,4) 
torus link. Similarly, the (10,4,3) resonance leads to two (5,2) torus knots, together forming a 
torus link of type (10,4). 

Note that the number n does not enter into the geometric specification of the knot or link. 



However, by (30), n mod k does determine the dynamics on the torus link. 






Figure 21: Schematic pictures of torus links of type (mi,m2) = (3,4), (4,4) and (10,4). Different colors 
distinguish the k — gcd(mi,m2) disjoint circles of the link. The central invariant circle (not part of the torus 
link) is also shown. 

Here we are specifically interested in local bifurcations at the elliptic invariant circle C = 



{r = 0} of (31 ). Some of these are analogous to the generic bifurcations of an elliptic fixed point 



of an area-preserving map in which periodic orbits are created or destroyed when the multiplier 
of the linearization passes through a root of unity. For the invariant circle, this would correspond 
to the transverse rotation number ujt passing through a rational, n/m2: an (0, rn-2, n) resonance. 
More generally, suppose that C is resonant, i.e., that a;(0, 0) satisfies a single resonance condition 



(26), and that 9^a;(0, 0) ^ 0. Then as /i passes through zero, a resonant torus will be created 



in the neighborhood of C. As above, this resonant torus is foliated into a one-parameter family 
of torus knots or links. If we now perturb (31), fe = /o + 0(e), then by analogy with the 



Poincare-Birkhoff theorem for area-preserving maps we expect that all but a finite number of 
these circles will be destroyed. Indeed, a generalization of this theorem to the volume-preserving 
map case has been obtained by Cheng and Sun |CS90b] . Thus the bifurcation should create a 
finite number of invariant circles in the neighborhood of C. 



Though the torus map (28) does not distinguish between its two angle variables, the embed- 
ding of (31) into M? assigns different roles to the longitudinal {9l) and transverse {6t) angles. 



The two-dimensional analogy then implies that for the volume-preserving case, ?tt,2 may play 
a more important role than mi. In particular, recall that an elliptic fixed point of an area- 
preserving map is "strongly" resonant if the denominator of ut is small, i.e., ra2 < 4. An 
analogous classification will pertain to the volume-preserving case. 

From these qualitative considerations and our numerical observations, we propose the fol- 
lowing conjecture: 

Conjecture 1. An elliptic invariant circle of a volume-preserving map with frequencies u) = 
{ujl,^t) that satisfy the single resonance condition coLmi +10x^2 = n (gcd(mi,m2,n) = 1, 
gcd{mi,m2) = k) generically undergoes one of the following bifurcations: 
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• 1712 = (String of pearls bifurcation): The circle is destroyed and a pair of saddle, period- 
mi orbits are born in a SCNS bifurcation, recall §5.^ The one- dimensional invariant 
manifolds of the saddles nearly coincide along the location of the destroyed circle. If the 
bifurcation is supercritical, it also creates a family of mi, almost-invariant balls (pearls) 
bounded by the manifolds of neighboring points on the orbits and containing elliptic circles. 

• |"T'2| = 1 (Saddle- center bifurcation): The transverse multiplier of the invariant circle 
becomes 1 and the circle is destroyed in a saddle-center bifurcation as the parameters cross 
the resonance line. 

• \i^2\ ^ 2 (Torus-link bifurcation): The invariant circle persists but may loose stability. In 
its neighborhood, k invariant circles are born that form an (mi, 7712) torus link. 

We believe that this classification of different types of m2-tupling bifurcations of invariant 
circles is new. In the following subsections, we will separately treat each of the three cases of 
the conjecture. 

7.1 Torus-Link Bifurcations, |m2| > 2 

In the beginning of this section, we discussed several examples of period doubling with \m2\ = 2, 



showing that a torus knot (link) is created if mi is odd (even), recall Fig. 20 We will now present 



more examples with larger \m2\, and then give an explanation of the observed bifurcations. 



Several of such bifurcations for /i = —2.4 are illustrated using Poincaree slices in Fig. 22 A 
slice is analogous to a Poincare section of a three-dimensional flow; however, since an orbit of a 
three-dimensional map would almost never intersect a two-dimensional section, we must instead 
consider a slice of some nonzero thickness. In the figure, we use the same slice, Sa = {{x, y, z) : 



y > 0, \z\ < A}, that was used to find the elliptic circles in ^6.1 

For example, at (e,/u) ~ (0.121,-2.4), the elliptic circle crosses the (4,-4,1) resonance. 
Since (mi, 771,2) have a common factor, this bifurcation creates an elliptic torus link of type 



(4,-4) which consists of four torus knots of type (1,-1), recall the schematic Fig. 21 The 
new circles move away from the central circle as e grows. In the Poincare slice — the top row 
of Fig. [22] — this bifurcation looks just like a generic quadrupling bifurcation of an elliptic fixed 
point in an area-preserving map. The slice reveals that both an elliptic and a hyperbolic (4, —4) 
torus link are created, and that the central circle appears to persist and remain stable through 
the bifurcation. We observe qualitatively similar Poincare slices for the (3,4,1) resonance (at 



(e,/i) ss (0.147,-2.4), recall Fig. [T6|), and (10,4,3) (at (0.169, -2.4)); however, the corresponding 
torus-links are quite different, recall Fig. |21[ 

It is easy to observe many such resonances; a list of the lower order resonances encountered for 
/i = —2.4 is given in Table|3]in Appendix [E) It is interesting to note that many resonances appear 
in sequences: (3, 4, 1) and (10, 4, 3) are the two first elements of the sequence (3 -|- 7n, 4, 1 -|- 2n). 
These resonances lines have one conjugate point, {loi,ujt) = (fjj^); iii common. Since this 
limit point is close to the image of the frequency map for fi = —2.4 many of these resonances 
are encountered in the family. However, since the slope of these resonance lines is proportional 
to n and the last observable, stable invariant circle has ujl ~ 0.2853 = [0, 3, 1, 1, 44, . . .], which 
is close to the rational y, the frequency curve will eventually miss the resonance lines. Any 
low-order, doubly- rational point in the image of the frequency map can similarly serve as a limit 
point for many resonance sequences. In the present case we observe, for example, the sequences 
(7n,l,2n), (7n, 2,2n), (3 -|- 7n, 3, 1 -|- 2n), and (4 -|- 7n, — 3, 1 -|- 2n). As a consequence, the 
invariant circle will repeatedly undergo bifurcations that have the same m2 but that correspond 
to different associated torus links. 
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Resonances with larger m2 have similar structure: for example a (1,-5,0) resonance at 
(e, /x) ~ (0.026, —1.1) creates a new elliptic and a new hyperbolic circle that, as e grows, give 
rise to a five-island structure when viewed in a Poincare slice. 

Tripling bifurcations {\m2\ = 3) may result in a brief loss of stability of the invariant circle. 
For example, there is a gap in the // = —2.4 frequency map of Fig. 16 near e = 0.235 where the 



circle crosses the (3,3, 1) resonance. This bifurcation again resembles the two-dimensional case. 
Indeed the new hyperbolic and elliptic (3, 3) torus links are created in a circle-saddle-center 
bifurcation before the "tripling", near e = 0.234. Moreover, the central circle appears to lose 
stability in a collision with the hyperbolic period-three family near e = 0.235; cross sections are 

3, 1) resonance occur for a range of /i, as can be seen 

0.252; this is a (10,3,3) 



shown in Fig. 22 This and the similar (3, 

-2.4 occurs in a gap around e 



in Fig. 19 Another tripling for // 



resonance. For this case {1711,1112) are coprime and the bifurcation results in the formation of 
a single, elliptic torus-knot with three longitudinal wraps and ten meridional wraps (the dark 
blue circle in Fig. \16[ . Similar tripling bifurcations occur whenever ^ > —1.6 as the invariant 
circle crosses the (1,-3,0) resonance. However, for these cases the tripled circle is a (1,-3) 
torus knot. 
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Figure 22: Poincare slices Ea for /i = —2.4 and the e values shown. The upper sequence crosses the (4, —4, 1) 
resonance, the middle the (3,3, 1), and the lower the (10,3,3) resonances. 

From the examples discussed so far it seems as if the bifurcations of invariant circles are 
simply those of fixed points of area-preserving maps multiplied by a (twisted) circle. However, 
this appears to be false. We found many portraits in the Poincare slice that are not generic in 
area-preserving maps. To discuss this, we first recall the generic bifurcations of a fixed point of 
an area-preserving map. When the multiplier of the fixed point is exp(ib27rin/m), n, m coprime, 
then there are 5 cases (assuming the generic higher order transversality conditions are satisfied) 
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|MH92] : 

• m = 1 (saddle-center): collision of saddle and center fixed points (double multiplier 1, not 
semisimple) , 

• m = 2 (doubling): a periodic orbit of period two is created and the fixed point changes 
stability (double multiplier —1, not semisimple). 

• m = 3 (tripling): a saddle-center bifurcation of creates a pair of period three orbits away 
from the fixed point, and the saddle passes through the fixed point, causing a momentary 
loss of stability. 

• m = 4 (quadrupling): this bifurcation can behave either like m = 3 (the hyperbolic orbit 
colliding with the fixed point with loss of stability) or like m > 4 (elliptic and hyperbolic 
pair being created without loss of stability) , depending on the relative strengths of higher 
order terms. 

• m > 4 (m-tupling): a pair of elliptic and hyperbolic periodic orbits of period m is created 
at the origin without loss of stability. These orbits form an island chain with m islands. 

As a result resonances are usually grouped into "low order" or "weak" (m < 4) and "high order" 
or "strong" (m > 4) cases. High-order resonances do not cause loss of stability of the fixed point 
at the bifurcation. 

For m > 3 these resonances are described by an area-preserving map of the form 

/ = e2™/"^z + 0(2) . 

The resonant terms correspond to the monomials z^ z "^ with ki,k2 G N and \ki — k2\ divisible 
by m. Upon setting z = y/zl exjpicj), for angle- action variables {(p,I), these bifurcations can be 
described by an approximate Hamiltonian 

H{(j), I) = fial + /3l'^ + ■ ■ ■ + {-f + 61 + . . .)I"'/2 cos mc/) , (32) 

where the bifurcation parameter /u detunes the resonance. When m = 3 the resonant term 7 
is the leading-order nonlinear term, and thus the behavior is different than for m, > 4. When 
m = A the leading-order resonant and nonresonant terms have the same order and the behavior 
depends on whether \j3\ ^ |7|. For m>b the resonant term is not of leading order and the usual 
island chain is born. 

Resonant bifurcations of invariant circles of volume-preserving maps follow a similar scheme, 
but there are more possibilities. This is for two reasons. Firstly, instead of the single period 
m there is now a vector (7711,7712). Roughly speaking, an area- preserving map describes a slice 
of the dynamics transverse to the invariant circle with "period" 7712. In the simplest situation, 
these 2Z?-dynamics rotate as the slice moves along the invariant circle; this is described by 7771. 
Secondly and less obviously, the distinction between low-order and high-order resonances is not 
so clear cut: indeed all the "low-order" resonances can also appear without loss of stability. Our 
observations can be summarized in three cases: 

• 7772 = 2: the central circle looses stability and an elliptic doubled circle appears. Depending 
of A; = gcd( 7771, 7772) there is either one circle that is a double cover of the original one 
{k = 1), or there are two circles each covering the original circle once that are mapped 
into each other {k = 2). These two circles form an (7771,7772) torus link. However, unlike 
generic bifurcations in Hamiltonian systems this bifurcation may be "weak", creating a 
pair of elliptic and hyperbolic circles. This happens, e.g., for the (7,2,2) resonance at 
(0.0167,-2.4). 
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• 1712 = 3: the tripling may lead to a momentary loss of stability of the central circle. When 
mi is divisible by 3, then 3 circles are created, and n mod 3 determines how they are 
mapped to each other. Unlike the generic area-preserving case, the tripling bifurcation 
can be "weak". This occurs, e.g., for the (11,-3,3) resonance at (0.1791,-2.4). 

• 1712 ^ 4: here a pair of elliptic/hyperbolic torus links of type (mi,m2) is created. If mi 
and 7712 are coprime, there is only a single circle, while if gcd(mi,?7i2) = k then there are 
k disjoint circles. This difference is invisible in a slice, and these bifurcations do look like 
generic bifurcations of area- preserving maps. The order in which the k circles are mapped 
into each other is determined by n mod k. 

It is possible that the additional cases when m2 < 4 could be described as non-generic 
bifurcations of area-preserving maps with m = m2 but in which the lowest order resonant term 



vanishes. The approximate Hamiltonian would still be given by (32), but now 7 = 0. For m = 2 
there are now two cases, depending on whether |/3| ^ \5\. Thus m = 2 becomes similar in some 
sense to m = 4 for the area-preserving case. When |/3| > \5\ an elliptic/hyperbolic pair is born 
like the generic case for 77i > 5. However, when |/3| < \6\, the bifurcation behaves like the generic 
m = 3 case: i.e., a hyperbolic orbit passes through the center with loss of stability. We have 
not yet observed this case, but we expect that is is also generic for bifurcations of invariant 
circles in volume-preserving maps. If 7 = for ttt, = 3, an elliptic/hyperbolic pair of orbits are 
born: this case behaves like the generic case m > 5. The mechanism by which the lowest order 
resonant terms vanishes in these cases presumably is discrete symmetry. We hope to return to 
this question in a future work. 

Besides the local bifurcations discussed so far there are also some interesting sequences 
of bifurcations associated with period doubling. We have already discussed several doubling 
bifurcations {\m2\ = 2) that result in the permanent loss of stability of the invariant circle, and 
thus indicate the loss of most of the bounded orbits. There are also doublings for which C is 



unstable only for a small range of e. For example the frequency map curves shown in Fig. 19 
sometimes pass through the (4, ±2, 1) and (5, ±2, 1) resonance lines — with a short interruption — 
indicating that the original circle regains stability after creating the doubled circle. Even when 
I mi I is large, a doubling bifurcation can result in the temporary loss of stability of C. For 



example, the gap in the frequency map of Fig. 16 near e = 0.302 corresponds to a (46, —2, 13) 



resonance. A pair of stable, circles is created that each wrap around the longitudinal direction 
once, while encircling the original circle in the transverse direction 23 times, an enlargement is 



shown in Fig. 23 



The class of torus knot (or torus link) bifurcations of invariant circles of volume-preserving 
maps of M^ shows a rich bifurcation structure. While there are similarities to area-preserving 
maps, there is also non-trivial three-dimensional structure and slices reveal some bifurcations 
that would be considered non-generic in the area preserving category. A more detailed investi- 
gation of these bifurcations would be very interesting. 

7.2 String of Pearls Bifurcations 

Resonances of an invariant circle with 7712 = correspond to rational values of the longitudinal 
rotation number. We will argue below that this bifurcation corresponds to a SCNS bifurcation 



for the map /™^, i.e., the bifurcation discussed in ^5.2 In the supercritical case the invariant 
circle typically degenerates into a pair of period-mi saddle-foci whose two-dimensional stable 
and unstable manifolds enclose a periodic family of elliptic invariant circles. The result is a 
structure that looks like a string of pearls on a necklace. This is a common phenomenon in our 
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Figure 23: Portions of the elliptic invariant circles near the (46, —2, 13) resonance when fi — —2.4. In the gap 
0.3017 < e < 0.3045 the central circle is unstable. The circles are shifted vertically for clarity. 



quadratic family, and we expect this to be a generic bifurcation of invariant circles of volume- 
preserving maps of M'^. Moreover, a similar bifurcation occurs in dissipative three-dimensional 
maps jBSVOSj ; they studied in detail an example with wl = I. 

For example when when e = and // = —2.0 for the map ([I|, w^ = 3;. As e is increased from 



zero, there is a period-four string of pearls created as shown in Fig. 24 In the figure the new 
period-four family of invariant two-tori and several orbits that spiral around the one-dimensional 
manifolds of the period-four saddle-foci are clearly visible. 



(e, [i) = (0.01, -1.9) 



(e, n) = (0.08, -2.0) 



"■" % 







V--' 



Figure 24: Phase space near the lol ~ \ string of Pearls bifurcation. The red cube is centered at the origin and 
has sides of length ^ y^. 

This structure can be expected generically. Suppose that as a;^ ^ ^ the invariant circle 
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C remains elliptic and that / can be approximately described by the SCNS normal form (18) 



When (jJl is rational every point on the invariant circle of the normal form becomes ?7ii-periodic 
with multipliers (1, A™^, A"'"^), where A = e^'^*'^^. In the resonant case, there are additional 
terms that must be kept in the normal form, but if |mi| > 3 and lot is irrational, then the 
quadratic terms in / do not change. Consequently /"^^ corresponds to the critical case on the 
saddle-center-Hopf line of Fig. [T] and, when any parameter is varied, it will generically undergo 
a SCNS bifurcation. If this bifurcation is supercritical, the singular fixed point of /"^^ will be 
replaced by a pair of saddle-focus fixed points and a new elliptic invariant circle will be created. 
For the original map this corresponds to the creation of a pair of period-mi orbits and a period- 
rrii family of invariant circles linked with the original circle. Recall from Appendix [D] that when 
the unfolding of this bifurcation is constructed, the coordinate transformations include a shift 
along the direction of the unit eigenvector: this shift selects the particular points along C that 
will correspond to the fixed points. 
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Figure 25: String of Pearls bifurcation for jj, — —1.383 and the values of e shown, 
vertically and scaled by £~2 so that their size appears to be roughly fixed 



The circles arc shifted 
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When e = 0, this 



Another example of this bifurcation is shown in Fig. 
longitudinal rotation number occurs at /x ss —1.3819 according to ([7]). Fixing ^ = —1.383 and 
increasing e from zero, the (5,0,1) resonance occurs at e ~ 0.069 where the "circle" develops 
five small bulges formed from the nearly-coincident, two-dimensional manifolds of two new, 
period-five saddle-foci. As e increases the points on the new periodic orbits move apart and the 
corresponding pearls grow in size. The one-dimensional manifolds of the saddle-foci also nearly 
coincide, forming an approximate circle up to e ~ 0.09. The string of pearls is still enclosed 
by a family of two-tori; one of these is shown in the uppermost portrait at e = 0.11. These 
tori and the new period-five tori apparently vanish near e = 0.12 and when e reaches 0.18 the 
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computation of q3reports no bounded orbits. A similar scenario occurs at (e, fi) = (0.076, —2.6) 
c„„ T?;„ Oft jjj ^]^jg gggg ii^Q pearls exist only for a very small window in e, and 



where ujl 



10 



see Fig. 
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the circle reforms as e is increased. Note that the ellipse algorithm is typically able to continue 
the circle through the string-of-pearls bifurcation when the slice is not too close to the periodic 
orbits. 




Figure 26: String of pearls bifurcations for /j, = —2.6 (lower three portraits) and /i — —3.5 (upper three 



portraits), where ul 



and I, respectively. 



In some cases an (?7ii, 0, n) resonance is responsible for the final destruction of the invariant 
circle. This happens for /i = —3.5 where the circle is undergoes an (8, 0, 3) bifurcation at 
£ = 0.266, see Fig. |27j A sequence of three portraits for this case is also shown in Fig. 26 
As for the g and -^ cases, an approximate, frequency-locked invariant circle appears to persist 
for a small parameter interval beyond the bifurcation (at least up to e = 0.2686); however, in 
this case the circle does not reconstitute as e increases. Destruction of the invariant circle also 
occurs for fi = —2.4 when w^ = ^ at e = 0.32. This bifurcation was shown in the frequency map 



Fig. 16 The (7, 0, 2) bifurcation appears to be subcritical as we can find no stable orbits beyond 
e = 0.32. This is suggested by the rightmost pane of Fig. [8] where the first iterates of seven 
unbounded orbits are plotted. These orbits appear to be near the two-dimensional unstable 
manifolds of a period-seven saddle and these manifolds do not appear to intersect as would be 
predicted by the subcritical case in Fig. [6| 

7.3 Pairs of Saddle-center bifurcations of invariant circles 

As we have seen, resonances on an elliptic invariant circle can lead to a local bifurcation that 
creates new invariant circles linked with the original. Invariant circles can also be created by 
resonances on two-tori and more generally by saddle-center bifurcations that are not associated 
with any regular structure. Sometimes saddle-center bifurcations are related to other resonances. 
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Figure 27: Frequency map for /i — —3.5. Inserts show the (8, —2, 3) resonance at e = 0.12 resulting in a doubled 
circle, (7, 12, 3) at £ = 0.181, (8, 2, 3) at £ = 0.22 and the final disintegration of the invariant circle in the (8, 0, 3) 
resonance at e = 0.266 leading to the formation of a string of pearls. 

For example, near the tripling resonance with m2 = ±3 (if 7 7^ 0) a pair of new invariant circles 
is born in a saddle-center bifurcation, as seen in Fig. [22| In the following we discuss a new 
saddle-center bifurcation that is caused by a resonance with 1712 = ±1. 

We observe that these circle-saddle-center bifurcations come in pairs, as can be seen in the 
frequency maps of Fig. 17, for 1712 = ±1. In these cases, the circle C created at e = approaches 
a resonance line of the form (?Tii,ibl,n) and its frequency map curves dramatically to "avoid" 
crossing the resonance. As this happens, an elliptic/hyperbolic pair of circles are born in a circle- 
saddle-center bifurcation far from the central circle. This may happen inside a chaotic region or 
region of unbounded motion or within the family of elliptic two-tori surrounding the invariant 
circle. In the latter case a torus bifurcates into a figure eight crossed with a circle. According to 
our observations, the rotation numbers of the elliptic member of the new pair lie on the opposite 
side of the resonance line. When ^ = —1.5, this is the (5, —1, 1) resonance line shown in Fig. 17 
A sequence of Poincare sections for this case is shown in Fig. [28] The new circles appear near 
£ = 0.097; they form a (5, —1) torus knot. Note, however, that the new elliptic circle (red in the 
figures) is relatively smooth and it is more appropriate to think of the highly deformed original 
circle (blue) as a (5, —1) torus knot wrapping about the new one. As e grows the stable region 
around the original circle shrinks and that around the newly created elliptic circle grows until 
it takes over as the primary elliptic circle and persists up to e = 0.175. The old orbit (blue) 
eventually dies in a circle-saddle-center bifurcation. So overall there is a pair of saddle-center 
bifurcations of invariant circles, first the creation of a new elliptic circle, and later the destruction 
of the old elliptic circle. This bifurcation also occurs at (e,/i) = (0.049,-1.1) corresponding to 
the (6, —1, 1) resonance. When fj, = —1.9 a similar bifurcation appears to occur near e = 0.07 
for the resonance (4, 1, 1), however, there are also many other resonances near these parameter 
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values that complicate this case. 
E = 0.090 
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Figure 28: Circle-saddle-center bifurcation for p = —1.5. Five Poincare slices show the creation of a new elliptic 
circle with (5, —1) helicity near e = 0.096. The stable region around the new circle grows while that around the 
original circle shrinks, until near e — 0.101 the original circle looses stability. Tori near the new original circle are 
shown in green and near the new circle are shown in red at e = 0.099 in the last inset. 



8 Conclusion 

The three-dimensional analogue of the two-dimensional, area-preserving Henon map can be 
found by normal form expansion and unfolding near a triple-one multiplier [DMO81 : we believe 
this quadratic map, M, should serve as the prototype for volume-preserving dynamics in M'^. 

We computed the volume of bounded orbits for ([I]) as a function of its two essential pa- 
rameters, (e, ^), obtaining an intricate structure reminiscent of Arnold tongues for circle maps. 
In the Henon map bounded motion appears near the elliptic fixed point. By contrast, in the 
quadratic volume-preserving map bounded motion appears near elliptic invariant circles. These 
invariant circles in turn are created in saddle-center-Neimark-Sacker bifurcations. These oc- 
cur on a CO dimension-one line in parameter space that emanates from the triple-one- multiplier, 
codimension-two point. The SONS bifurcation also creates two fixed points that are at the poles 
of the vortex bubble containing the bounded orbits. 

Near the SONS bifurcation, the map can be decoupled in cylindrical coordinates into a 
semi-direct product of an area-preserving map independent of the angle and a map of the angle 



depending on all coordinates, (18). A consequence is that close to the bifurcation, the two fixed 



points have one-dimensional manifolds that nearly coincide along the bubble's polar axis and 
two-dimensional manifolds that form the bubble's outer boundary. The curves of transversal 
intersections of the latter manifolds undergo bifurcations that were described in |LM98| ; these 
give rise to bounded chaotic motion near the boundary of the vortex bubble. 

When the SCNC bifurcation is supercritical, it creates an invariant circle C in the vortex 
bubble that corresponds to the third, elliptic, fixed point of the area-preserving, decoupled map. 
We computed the location and frequencies of C for the normal form and obtained analytical 
expansions for the longitudinal, ujl, and transverse, ujt, frequencies of the invariant circle (or in 
fact for any torus nearby) for the map M. 

The frequency map, 0, is an essential object of KAM theory. In a Hamiltonian setting 0, 
maps n actions into n frequencies for tori of maximal dimension. In the volume-preserving case 
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the dimension of the action and angle spaces are typicahy not equal. In M^, the tori in the neigh- 
borhood of the circle C have effectively one action, and Q defines a curve in the two-frequency 
space that depends upon that single action. As is done in the study of lower-dimensional tori in 
Hamiltonian systems, the frequency map can be extended to a local diffeomorphism by including 
parameters in the domain. In this way we study the frequency map 0, : (e,^) -^ (ujl,ujt) of 
the invariant circle of the quadratic map. This mapping organises the bifurcations of invariant 
circles. Our analytic approximation for Q compares very well with the numerical computations. 

In order to numerically compute O, we developed a simple method to find the location of 
elliptic invariant circles of maps. The idea is to consider a slice of the mapping that is transverse 
to the invariant circle and fit a thin ellipsoid to iterated points returning into the slice. This 
method converges very fast and appears to be a robust and accurate way to continue the invariant 
circles. Having found a good approximation of C, the computation of the longitudinal rotation 
number can be made quadratically convergent by considering its continued fraction expansion. 
For the transverse rotation number we have a linearly convergent numerical method. 

Equipped with the continuation method and the frequency map an astounding variety of 
resonant bifurcations of invariant circles can be found whenever rriiuji + m2UJT = n. We clas- 
sified these into three cases: string-of-pearls bifurcations (m2 = 0), paired circle-saddle-center 
bifurcations (m2 = 1), and torus link bifurcations (m2 > 2). More precisely a torus knot of type 
(m-i, m-2) is found when k = gcd(mi, 7712) = 1, while a torus link of type {mi, 1712) consisting of k 
circles is found when A; > 1. For example, there are two types of circle-doublings: either the new 
circle is a double cover of the original one, or there are two new invariant circles, each covering 
the original circle once. 

The most spectacular bifurcation is the string of pearls that occurs when the longitudinal 
frequency becomes rational. Here the invariant circle breaks into a pair of saddle-focus periodic 
orbits, and each pair recapitulates the original vortex bubble on a smaller scale. In this way, 
we expect that is is possible to find self-similarity in the SONS bifurcation: the invariant circle 
bifurcates to create new bubbles that contain invariant circles linked to the original one and that 
will themselves undergo string-of-pearls bifurcations under parameter variation, etc. It seems 
likely that the there exist parameter values with infinite sequence of linked pearls around pearls 
around pearls .... 

Similar bifurcations also occur in dissipative 3D maps. For example circle doublings can 
occur, but unlike the infinite sequences of doublings of periodic orbits, these appear to be 
typically finite in one-parameter families |AA84j . Moreover, a bifurcation similar to the string 
of pearls also occurs in the dissipative case near the codimension-two a saddle-center point. 
The three-parameter unfolding of this bifurcation is discussed by Broer et al |BSV08j ; they also 
studied a model of the dynamics near w/, = g. It is interesting that they observed that there 
can be an attracting invariant two torus containing a repelling period-five bubble. We hope to 
apply similar methods to study the volume-preserving, resonant cases in a future paper. These 
are both simpler (since there is one fewer parameter) and more complex (since there are no 
attractors) than the dissipative case. 

A Appendix: Parameter Reduction 

Though the map (IM nominally depends upon the three unfolding parameters {e, fii, ^2)^ one of 
these can generically be eliminated by an affine shift in x. If we set x ^ xq + x then the new 
map has the same form ^ with p{x, y, e, jjl) replaced by p{xq + x, e, fi). If there is an xq so that 
dxP{xo,0,£, fj.) = 0, then the linear term fiix can be eliminated. Near the origin, this equation 
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has a solution providing d^p{0, 0; 0, 0) / OJ^ In this paper, we assume this is the case, so that 
without loss of generality we can set fii = 0. Moreover, even though p is nominally a function of 
(e, fi) we will take the coefficients of this function to be independent parameters thus suppress 
this dependence. The normal form (|5| then becomes M. 

For the quadratic case M, one of the parameters in P could be eliminated by scaling the 
variables; for example, under the scaling ^ -^ y^|a|^, the parameter a is effectively replaced by 
±1. Thus for this case (II| reduces to a four parameter family. 

B Appendix: Basic Properties 

In this section we study some of the simple dynamical properties of Q, especially the fixed 
points and low period orbits. At first we let P be an arbitrary, smooth function that is 0(2) in 
the variables (x,y), and such that 

a=^P,,{0,0)>0, (33) 

which we used to eliminate the parameter /ii, see Appendix [A} 

When P is smooth, the map ([I]) is a diffeomorphism; indeed, its inverse can be easily com- 
puted 

(X -y + z \ 

y-z . (34) 

z + e- fi{y- z) - P{x -y + z,y-z)J 

The fixed points of M have the form (x*, 0, 0), where 

P(x*,0) = e . 



By (33), this equation has two solutions. 



Xl ^±^'-{1 + 0(8)) , 

that emerge from the origin for e > 0. This corresponds to the saddle-center bifurcation in the 
normal form. Indeed, the Jacobian of / at these points has characteristic polynomial (|3| with 

r^ = 3 + fi + Py{xl,0) , a± = 3 + ^l + Py{xl,0)-P^{xl,0) . 

The saddle-center line, corresponding to a multiplier Ai = 1, occurs on the line r = o" in Fig. [T} 
Since a ^ 0, this corresponds to the line e = 0. 

Recall that a diffeomorphism is reversible if it is conjugate to its inverse, r o f = f^^ o r. 
The normal form M is only reversible if P satisfies special conditions. For example, if 

P{x,y) = P{-x-y,y) , 

then there is an involutory reversor 

r{x,y,z) = {-x,y- z,-z) . (35) 



^ If this is not the case, but dyp{0, 0, 0, 0) 7^ 0, then we can choose xq so that dyp{xo, 0, e, fj,) =0 ehminating 
the hnear term in y. 
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This condition on P is essentially necessary for reversibility as well. Recall that a reversor maps 
each period-n orbit into another such orbit with reciprocal multipliers. Orbits that lie on the 
fixed set of r, in this case the y-axis, are mapped into themselves by r and are called symmetric. 
Orbits that are related by the reversor have reciprocal multipliers, and a symmetric orbit must 
have self-reciprocal multipliers and therefore one unit multiplier. When e 7^ 0, the fixed points 
do not lie on the saddle-center line and consequently cannot be symmetric. If there are only two 
such points, then they must be mapped into each other by the reversor and must have reciprocal 
multipliers. This means that t± = cr^:, or equivalently 

Py{xX,o) - Py{xr_,o) = p,.(x;,o) = -p^{x*_,o) . 

For the quadratic case, this implies that P has the form 

P = a(x^ -|- xy) + cy"^ . 

and that the reversor is (35). Consequently, the map (IT]) is only rarely reversible. 

For the case that P is a positive definite, quadratic form, that is o > and D = ac—b'^ /A > 0, 
it was shown in |LM98j that every bounded orbit of the quadratic map is contained in the cube 
l^loo < i^B where 

KB = a(\3 + fi\ + l + V(|3-F/i| + iy + \e\/a_ 



where a = jj max(c, a — b + c)\^^\ 

C Quadratic Map Fixed Points 

The stability properties of the fixed points (pi) of the quadratic map(IT| are easily determined 
from the characteristic polynomial (p|. For the fixed points the trace and second trace become 

T± = 3 + fx + bx± , 

a± = 3 + fj, + {b — 2a)x± . 

Since a > 0, the stability diagram Fig. II] maps diffeomorphically onto (/U, x-t) plane, the half- 
planes r — o" ^ corresponding to x±, respectively. As for the case of a general nonlinear 
function P, the saddle-center bifurcation curve corresponds to e = 0. Along this curve, when 
—1 < T = a < 3 the other two multipliers have the form e '^'^^^ which implies that —4 < /x < 
and Lo = ujq{^), (|7|. Thus we call this the saddle- center- H op f segment. 

Note that for most parameter values, one fixed point is type (2, 1) (two-dimensional stable 
manifold) and the other is type (1,2) (two-dimensional unstable manifold). The exception 
corresponds to a parabolic region near // = —4 that happens to coincide with the existence of a 
period-two orbit. Indeed, the period-doubling bifurcation in Fig. [T] corresponds to r -|- a = —2. 
For the quadratic map, this implies ^ = (a — b)x± — 4, which is equivalent to the parabola 

A = (a-5)2e-a(/i-h4)2 = , (36) 

as shown by the dashed lines in Fig. [2] It is easy to see that the quadratic map has at most one 
period-two orbit, which has the form (x, y, z) — > (x -|- y, — y, —z) -^ (x, y, z) where 

(^' ^' ^) = ^^3^ (m + 4 T /o, ±2p, ±4p) . 
^^ This corrects an error in |LM98| in the formula on page 572 for k. 
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Consequently, this orbit exists only when a^h and 

2 A 

For the case that a — 26 + 4c > 0, the orbit exists in the interior of the period-doubling parabola 
(the dashed curves in Fig. ^; otherwise, it exists in the exterior of this parabola. The charac- 
teristic polynomial for the period two orbit has the form (|3| with 

4A 
T — a = . 



As expected, the period-two orbit is born on the parabola A = with a multiplier 1. Thus, 
for fixed (a, 6, c) only four of the eight stability regimes of Fig. II] are reached; an example for 
a > b and a -|- 4c > 26 is shown in Fig. [29l In this case, the period-two orbit is created on the 
upper half of the doubling parabola A = by being emitted from the fixed point (x-|_, 0, 0) and 
is destroyed on the lower half of the parabola by being absorbed into the second fixed point 
(x_,0,0). 

The two curves of double eigenvalues in Fig. Ill cor respond to Ai = A2 = r giving r = 2r-|-l/r^ 
and a = r'^ + 2 jr. For the fixed point these become 

2ax± = -r^ +2r- 2/r + l/r^ , 

ji = hx± - ( 2r -I- ^ j - 3 . 

These are shown as solid curves in Fig. [2j 

Finally, fixed points with multipliers of the form A = [re^'"'''^ , re"^'^*'^, -j) have r = 2r cos(27ra;)-|- 
1/r^ and a = r'^ +2 cos(27rL<j)/r, which implies 

2ax± = ir I I 2 cos(27ru;) — r j , 

\ rj \ rj ^3^^ 

/i = 2r cos(27ra;) -\ — 2 ~ "^ ~ ^^^ • 

Several curves of this type for fixed, rational values of u; are shown in Fig. [2] 

The quadratic map ([I]) has at most one pair of period-three orbits; this follows because the 
set of equations ^ = /^(^ can be reduced to a degree-eight polynomial in one variable. This 
polynomial has a quadratic factor with zeros at the fixed points; the remaining sixth degree 
equation has zeros on the pair of period three orbits. They are created in a saddle-center 
bifurcation along the curve 



where 
K = 



a{2a -5b + 6c)((a - b)^ + 3c(a - b + c)) 
(a2 - 3a6 + 3ac + Sft^ - 96c + 9c'^){a^ - 3a% + 2,a?c + ac^ + 2b'^c - babe + 2,ab'^ - P) 



For example, when a = l,6 = c= 21-?^ = 4. Note that this curve is not the same as the tripling 
curve of the fixed point. 

One of the period-three orbits is type-(2,l) and the other type-(l,2). When /i < —3, 
the period-three orbits are born with real eigenvalues, corresponding to r = o" > 3, and (for 
a = 1, 6 = c = 2) when —3 < fi < —2.24 they are spiral-saddles (saddle-center-Hopf bifurcation) . 
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Figure 29: Stability diagram for the period-two orbit of (IIJ with the polynomial S and a = 1, b — c = 0.5. The 
dashed (black) parabola is the curve A = where one multiplier of the period- two orbit is 1. This is the same as 
the period-doubling curve in Fig. [2] The period-two orbit exists for A > 0, to the right of the A = parabola. 
The solid (green) curves correspond to a double eigenvalue, and the dotted (red) curve to a multiplier —1. 

At /i ss —2.24 the period three orbits sit at the point t = a = —1, the (—1, —1, 1) point, e just 
to the right of this parabola, these orbits have real eigenvalues. They then cross a curves of 
double multipliers, and become spiral saddles. 

Overall it should be noted that the importance of fixed pointed for three-dimensional, volume- 
preserving maps is far less than for two-dimensional, area-preserving maps, since they are gener- 
ically unstable. The important object in the 3D volume-preserving case are elliptic invariant 
circles. 



D SCNS Normal Form 

In this section we summarize the computation of the normal form for the saddle-center-Neimark- 



Sacker (SCNS) bifurcation of a volume-preserving map that is discussed in ^5.2 and thus the 
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technical part of the proof of Th. [1} In the complex coordinates (n, v, z) the map takes the form 

C' = /(C) = M{C + b{C)) , M = diag(e2--, e'^--, 1) , (38) 

where b = 0{2) represents the nonlinear terms. 

The formal, normal form for / is obtained by constructing a near identity transformation 
■q = '(/'(C) = C + ^(C) such that the conjugacy, g o ip = ip o f gives a new map 

g{7]) = M{7] + cir])) 

that is "simpler" than / — for example contains fewer nonlinear terms. Assuming that / is in 
normal form through terms of degree d — 1, this means that we must find a degree-d vector 
function h such that the conjugacy holds to degree d. To do this we must solve the homological 
equation, 'LMh{C) = c(C) — ^(C)i at degree d for h where 

LmHC) = M-'h{MC) - h{C) . (39) 

The terms in b that are in the range of Lm can be eliminated by the transformation. Thus the 
nonlinear terms c in the normal form correspond to the parts of 6 in a complement to the range 
of Lm- 

The normal form is formal in the sense that the map and transformation are expanded in 
formal power series in the variables {u,v,z). A basis for the space of vectors of polynomials of 
degree-d are the "vector monomials," 



P 



t,m 



Cci , niGN^ , (40) 



where Cj ,i = 1,2,3 are the unit vectors Cj , = 6ij, C™ = vJ^'^v'^'^z'^'^ and \m\ = S? iTtij 



d. When M is diagonal, the homological operator (39) has a diagonal representation in the 



monomial basis (40): 

^ - 1 j Pi,m ■ 

Consequently ker(LAf ) is spanned by those monomials that satisfy the "resonance conditions" 

A*" = A, , (41) 

and we can choose the normal form to lie in the kernel of La/. The elements of the kernel 



correspond to eigenvectors with zero eigenvalue; if we assume that w in ( 16 ) is irrational, then 
the resonance condition is satisfied for arbitrary integers m^ and when 







1, 


i = 1 


mi - 


- 77l2 = " 


-1, 


i = 2 






[ 0, 


i = 3 



Thus to arbitrary degree, the normal form can be written 

u = \u{\ + Cl(|ti| , z)) , 
z' = Z + C3{\u\'^ , z) , 



(42) 



where Cj are arbitrary polynomials subject to the requirements that Cj(0, 0) = and DzC^^O, 0) 
0. 
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Since / is volume-preserving, h can be selected to be volume-preserving as well |DM08] : to 
implement this efficiently, it is better to use Lie series instead of power series. We compute the 



power series by truncating the exponentials. Through cubic degree (42) becomes 

An(l + Az + Bz^ + C|u|2) + 0(4) , 



u 



z' = z + az"^ + /3|u|^ + 7|Mp2; + kz^ + 0(4) , 
where A,B,CgC and a, /3, 7, k G M. This map is volume-preserving to 0(3) when 



(43) 



-a 



2Br = 3(a^ - k) - 
4Cr = -7 - 2a/3 



A 



(44) 



where we denote A = Ay + lAi, etc. 
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The normal form (42) can be unfolded by assuming that the map / depends upon a set of 
parameters p. Following Elphick et al |ETB^87] , formally expand /, h, and g in double power 
series in p and ("; then, for each degree in p, the coefficients of h must still satisfy the homological 



equation. To zeroth order in p and all orders in C,, the normal form is given by (42). The new 



feature occurs at first order in p: there are now terms in / that are constant and terms that are 
linear in the variables C,. For the constant terms, the homological equation reduces to 



(M-/)/io = co-6o 



(45) 



Thus any terms in b not in the null space of M — / can be removed. This implies that there is 
generically an affine term (0, 0, 5{p)) in the normal form at this order. Since this is the sole term 
that occurs at this order, we can use (5 as a primary parameter; it corresponds to unfolding the 
saddle-center bifurcation. 

The terms of first degree in C, can be represented by matrices, e.g. /ii(C) = HC^ ^tc, and the 
homological equation becomes the matrix equation 



M-^HM- H = C -B 



(46) 



Since M is diagonal, this operator is diagonal in the components of H and its kernel consists of 
the matrices that commute with M: the diagonal matrices. Thus the normal form must include 
an added diagonal, linear perturbation. For the volume-preserving case when w 7^ or 
means that 



^, this 



M^diag(Ae^™,Ae 



-2mu) 



'A2' 



for a new frequency uj = uj + 0{p) and expansion/contraction factor A = 1 -|- 0{p). The three 
parameters (5, A,u;) unfold the saddle-center Hopf-bifur cation; however, as we will see below one 
of them can generically be eliminated. 

For all terms of higher degree in C,^ the homological equations are the same as before; con- 



sequently, the form of the terms in the normal form is the same as in (43), except that the 
coefficients are now allowed to depend upon the parameters. 
To quadratic degree, the unfolded normal form is thus 



u' = Ae^^'^u{l + Az) + 0(3) , 

z' = -5 + zjk^ + az^ + /?|u|2 + 0(3) . 

Alternatively, we could eliminate a — —Ar, 7 = — 4Cr + 2Arl3, and k = a^ 



(47) 



iBr 



1 A-^ 
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/" 


e 


mi 


m2 


n 


-1.95 


0.23 


4 





1 


-2.4 


0.318 


7 





2 


-2.6 


0.0760 


10 





3 


-2.6 


0.095 


283 





85 


-2.7 


0.027 


13 





4 


-2.8 


0.0354 


41 





13 


-3.3 


0.1164 


14 





5 


-3.5 


0.265 


8 





3 


-3.8 


0.24 


12 





5 


-3.85 


0.152 


7 





3 


-3.85 


0.305 


19 





8 


-3.9 


0.053 


9 





4 



Table 2: Parameters for some of the string-of-pearls bifurcations 



Note that when a is nonzero, the imphcit function theorem implies the existence of an affine 
shift z = z + /S. such that the new map has unit multiplier at the new origin. Hence this shift 



can be used eliminate A, and the final normal form is identical to (43) with the addition of an 



affine term —5 in the z' equation. This coordinate transformation modifies the rotation number 
u) as well. 



Transforming (43) to symplectic, cylindrical coordinates gives the final result (18). 



E Appendix: Resonance tables 

Table [2] shows some parameter values for different string-of-pearls bifurcations. Table |3] shows 
lower resonances encountered in the family with ^ = —2.4. 
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Table 3: Approximate resonances, \m- oj — n\ < 0.001, with |mj| < 10 for invariant circles with 
/i = —2.4, see table [T| page 20 and Fig. 16 
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